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Lower Case Latin 


a semi-diameter perpendicular to Z-axis in ellipsoid 
b semi-diameter parallel to Z-axis in ellipsoid 
f distributed loads in local coordinates 
h shell thickness; face sheet thickness 

i index: beginning edge of shell segment; independent joint of 

kinematic link; subscript "inside" 

j index: ending edge of shell segment; dependent joint of 

kinematic link 

n index on harmonic 

o subscript "outside" 

r radius 

s index of segment; coordinate in cylinder or cone 

t core thickness in sandwich shell 

w normal deflection, positive inward 

Upper Case Latin 


C stiffness eccentricity parameters; offset distance in ogive, ellipse 
D bending stiffness parameters 

E Young’s modulus (lb/in ) 

F lineal force (ib/in) 

G shear modulus (lb/in ) 

K extensional stiffness parameters 

M bending moment on shell (in-lb/in) 

N membrane stress resultant (lb/in) 

N assumed membrane stress resultant (lb/in) 

Q transverse shear stress resultant (ib/in) 

R radius; "global" coordinate, positive radially outward. 
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temperature; "global" coordinate, tangential 
Cartesian coordinate, 6 = 0 at X-axis 
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Cartesian and. "global" coordinate, coincides with axis of revolution 

Greek 

angle between rotated coordinates 

ratio of semi-diameter parallel to Z-axis in ellipsoid to semi- 
diameter perpendicular to Z-axis 

shear strain; non-linear parameter; .angle of inclination of 
kinematic link 

normal coordinate, positive inward 

circumferential angular coordinate (rad) ' 

shell parameter 

Poisson’s ratio 

normal stress (lb/in ) 

shear stress (lb/in ) 

meridional angular coordinate (rad) 

rotational displacement (rad) 

rotational displacement in "global" coordinates (rad) 
displacements in fixed or "global" coordinates 
segment length parameter 

Miscellaneous 


eq equivalent 
s<f> sin 
Ctj) cos <j> 

Other symbols are defined in the text where used. 
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SECTION 1 


PROGRAM CAPABILITY 

< 

The use of an accurate shell theory to analyze structural shell problems 
usually involves complex mathematics and numerical techniques, which are 
nearly impossible to treat without the aid of automated procedures. On this 
basis, a digital computer program based upon the Love-Reissner first order 
shell theory has been developed. This program can analyze orthotropic thin 
shells of revolution, subjected to unsymmetric distributed loading or con- 
centrated line loads, as well as thermal strains (Reference l). Furthermore, 
a shell with arbitrary boundary conditions, under loads which vary arbitrarily 
with position and under a temperature variation through the thickness, is 
tractable with this program. The shell can consist of any combination of the 
following geometric shapes: 

1) Ellipsoidal - spherical (offset from the axis of revolution allowed) 

2) Ogival - toroidal 

3) Modified ellipse shape 

4) Conical - circular plate 

5) Cylindrical 

6) General point input geometry 

7) Dummy geometry slot to be filled by the user * 

8) Discrete ring 

9) Elastic support 

The shell wall crossection can be a sheet, sandwich, or reinforced sheet or 
sandwich. The reinforcement can consist of rings and/or stringers, a waffle 
construction rotated at any angle to the principal coordinates, or an isogrid 
construction. General stiffness input options are also available. The rein- 
forcement material properties can differ from those of the main shell, and a 
temperature variation can cause different properties in the two face sheets 
of a sandwich shell. 
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The basic approach to the problem (Reference l) is to cut the structure into 
several shell regions. These regions need to be singly-connected shells, 
and can only have line loads applied at their end points. There are no 
restrictions on geometry, or uniform or thermal loads. The regions are 
further subdivided into several shell segments, each being free to have its 
own geometric shape, provided that the shape falls into one of the categories 
mentioned above. 

Stiffness matrices obtained for each segment, are coupled by standard matrix 
methods to obtain region stiffnesses, which, after being reduced in size, are 
in turn coupled to form the total shell structure under analysis. Currently, 

the UNIVAC 1108 computer program is sized to handle a structure composed of 

< 

up to 2k- segments in each of 19 regions arbitrarily connected to each other. 
There is a limitation on the size of a shell segment, which is a consequence 
of the demand that boundary disturbances be felt throughout the segment. This 
limitation is mathematically described in Section 2 (pages 2- 31 to 2- 33) as a 
length parameter. This parameter, however, is not reliable near the apex of 
any shell shape ( <f> = 0), and the segments needed in this region are actually 
much smaller than predicted by the parameter. A mathematical, singularity 
occurs at the apex where t q (the radius of revolution) becomes zero. It is 
this singularity which prevents the length parameter from being meaningful 
near the apex. Furthermore, the point (<j> = 0) is not an acceptable input 
point of the program (except for the torus-ogive and offset ellipsoid), 
although any, point outside a circle of infinitesimal radius is satisfactory. 

There is a considerable latitude in what can be done within each shell seg- 
ment. The thickness of any segment can be symmetrically tapered and it can 
contain up to l4 points of discontinuity, provided that the segment center- 
line remains continuous and describable by a single shell geometry. A temp- 
erature distribution through the thickness can be specified at three points 
in a homogeneous shell, and 4 points in a shell of rigid core sandwich con- 
struction. The distribution is considered to be linear between these points. 
Thus, it is possible to approximate temperature distributions other than 
linear distributions. In the event of physically discontinuous shell center- 
lines, a kinematic link is available for use in the analysis. The link 
relates displacements across the discontinuity. This link may be used between 


1-2 



regions, and between segments within a region. Discrete offset rings are 
also available for use within or between regions. 

The present program is also capable of a non-linear analysis of axi symmetric 
ally loaded shells. The analysis of this large deformation case is accom- 
plished by use of iteration. Details of the non-linear theory involved are 
presented in References 1 and 2, and the program utilization of this option 

4 

is described in detail in Section 2 of this report. In the linear analysis 
of unsymmetrically loaded shells, the partial differential equations of; the 
shell theory are reduced to ordinary first order differential equations, 
which are solved with the aid of a Runge-Kutta method of numerical integra- 
tion. The reduction is accomplished by use of a Fourier series expansion in 
the circumferential coordinate, 0. Axisymmetric loading is represented by 
the "zeroth" term of the Fourier series alone. For distributed loads, such 
as aerodynamic pressure, inertia loads, and aerodynamic heating, the Fourier 
series expansion is convenient since these loads vary smoothly in the circum 
ferential direction. For the most general inertia loads, only the zeroth, 
first, and second terms of the series are necessary. For reasonably smooth 
pressure distributions, the same terms will usually describe the loading 
adequately for the purposes of the stress analyst. Concentrated point loads 
can also be described by Fourier series expansions. Although, in general, 
many terms are necessary for a good representation, the SATELLITE - 2S pro- 
gram (Reference 3) is available for automated computation of these terms. 

Fundament ally different types of loading are represented by the various 
"harmonics". Pressure and temperature affect all harmonics. . In ’addition, 

i 

certain components of the three harmonics n = 0, 1, 2 reflect particular 
physical loads (see Reference k , Appendix A). v 

• n = 0 (Zeroth Harmonic): This is the axisymmetric contribution. 

Loads are due to axial translational acceleration, and centrifugal 
forces. Net axial load is produced only by this type of loading. 

• n = 1 (First Harmonic): This is the antisymmetric contribution. 

Loads are due to angular, centrifugal, and lateral translational 
accelerations. Net lateral load and a bending moment, which do not 
cause circular cross-sections to deform, are produced by this load 
contribution. 
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• n = 2 (Second Harmonic): Because of the character of the deforma- 

tion, this is called the "ovalizing" contribution. Loads are due 
to centrifugal accelerations (see Reference 4, Appendix A). 

No net loads can be obtained for n > 2; these higher harmonics contribute 
self-equilibrating loads due to pressure, temperature, and concentrated 
forces. For unsymmetric loading, the program will apply and sum all desired 
Fourier harmonic terms. A value V at a particular point (r Q> 0) on the shell 
is thus given by: 

N 

V(r 0 , 0) = ( r Q ) cos n9 

n = 0 

The user of the program has the option of obtaining only the final summed 
answers V(r , 0) at specified values of 0, or intermediate output from each 
harmonic also. The output is printed out for each segment of the shell at 
intervals specified by the user. In addition, an option for graphical dis- 
play of each segment stresses and displacements is available. 

In the case of a premature program termination due to time estimates in the 
execution of an unsymmetric loading problem, a restart capability is provided 
so that previous computations will not be wasted. 


; 
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SECTION 2 


INPUT INFORMATION 


2 . 1 GENERAL NOTES 


The preceding section provides some insight into the capability of the pro- 
gram, and the potential that it might have for future use. If the program 
is applied judiciously it can he an extremely powerful tool. The mechanics 
of applying it should be clearly understood. With this in mind, the re- 
maining section should be studied carefully. 

The required input data may be subdivided into three main parts, namely: 
geometric, topological (or coupling orientation) and joint data (degree of 
freedom description for each joint component). Each segment requires its 
own geometric configuration and numerical integration control. 

The output consists of stiffness coefficients for each shell segment and the 
actual symmetry of the coefficients is presented in a convenient form for a 
check on the accuracy of the integration through the segment. Region stiff- 
nesses and their symmetry checks are also provided. Final stresses, dis- 
placements and Huber-Von Mises-Hencky "effective stresses" are printed out 
for each shell segment at intervals along the segment as specified by the 
user of the program. The output will be further discussed in Section 3- 

The present program size is described in the table below. 

UNIVAC 1108 Table of Program Sizing 


I. 

Segments per region: 

2b 

II. 

Segment joints per region: 

25 

III. 

Regions: 

19 

IV. 

Region joints: 

20 
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UNIVAC 1108 Table of Program Sizing (continued) 


V. Number of points available 
per segment for specifying 
geometric or load data: 30 

VI. Number of points available 
through the thickness for 
specifying temperature data: 4- 

VII. Geometries: ellipsoid, sphere, offset ellipsoid, 

modified ellipsoid, ogive, toroid, 
cone, annular plate, cylinder, general 
geometry, ring, elastic support, dummy 
geometry. 

VIII. Wall cross-section options: single sheet, equal face sheet sand- 

wich, unequal face sheet sandwich, 
eccentric reinforcement (rings, 
stringers or both), waffle reinforce- 
ment rotated at an arbitrary angle to 
coordinate axes, isogrid reinforcement, 
arbitrary stiffness input. 

IX. Number of material property 

tables per submission: 10 

•X. Number of points per 

material property table: 10 

XI. Number of consecutive load 

conditions per submission: 10 (Except when summing harmonics or 

non-linear^ = l) 

XII. Orthotropy options: isotropic or orthotropic sheet, iso- 
tropic or orthotropic sandwich, iso- 

tropic or orthotropic sandwich with 
different face sheet properties caused 
by thermal gradients, isotropic or 
orthotropic sheet or sandwich reinforced 
by different property rings or different 
property stringers or both, isotropic or 
orthotropic sheet or sandwich reinforced 
by a different property waffle system 
rotated by an angle P to coordinate axes, 
isotropic or orthotropic sheet or sand- 
wich reinforced by a different property 
general isogrid reinforcing system, 
arbitrary stiffness input options used to 
describe other configurations. 

XIII. Number of harmonics in a 

multi -harmonic submission: 25 
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UNIVAC 1108 Table of Program Sizing (continued) 


XIV. 

Number of Q angles for 
description of circum- 
ferential variation: 

36 

XV. 

Number of rings per region: 

23 

XVI. 

Total number of rings at 
region joints in structure: 

18 

XVII. 

Number of elastic supports: 

20 

XVIII. 

Number of graphical display 
items : 

56 


Figure 2-1 shows the detailed option flow chart for the present program. 


GENERAL NOTES - Idealizations 


Before discussing the specific card input order, it would be advantageous to 
introduce some general guidelines in the area of idealizations and topology. 

In many computer programs there is such an abundance of numerical computation, 
that minimizing numerical roundoff errors becomes as important as getting the 
final answers. In some cases the engineer can aid the program in this effort 
through the use of judicious idealizations. Such a possibility exists in the 
STARS -2 programs, since many internal operations are involved with building 
and inverting stiffness matrices. The object of the user therefore, should 
be to help the computer by avoiding the creation of ill-conditioned matrices 
at any step (see Reference 5)* Physically, the way to achieve this end is to 
have all the segment stiffness matrices of the same order of magnitude. This 
will in turn produce region stiffness matrices which are of similar orders of 
magnitude, and minimize possible ill-conditioning in the total structure 
matrices. The user can help to achieve this end by sizing his segments in 
such a way so that no short stiff segment is contained alone in a region with 
all other long flexible segments, or that no region comprised of all short 
stiff segments exists in a structure whose other regions contain only long 
flexible segments. No accurate measure can be given on the relative stiff- 
ness or flexibility of segments allowed, and thus the best check is to see if 
a structure is in equilibrium under the applied loading. It also must be 
kept in mind, that if an idealization has provided useful results for 
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Figure 2-1. Program Option Flow Chart 
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axisymmetric (zeroth harmonic) loading, it need not necessarily be a good 
idealization for other harmonics, but if good results are obtained for the 
first harmonic, then other harmonics should also encounter no adverse 
behavior. The symmetry checks of segment and region stiffness matrices are 
useful for many reasons, but will not necessarily alert a user to ill- 
conditioning. 

In the use of regions, one other type of accident must be avoided. This is 
the creation of a single region structure with both ends fixed, wherein no 
suitable boundary condition matrix can be formed. Thus, in the use of region 
idealizations, which are less physically meaningful to a user than pure seg- 
ment idealizations, care should be taken so that al 1 boundary conditions are 
not zeroed out. To avoid this problem, and to minimize program running time, 
it is best to maximize the number of regions in a structure, and minimize 
the number of segments per region. Thus, in small problems, for best 
numerical efficiency, there should only be one segment per region. 

In the solution of problems involving only axisymmetric loads (n = 0) and no 
torsion, the torsional degrees of freedom should be removed. This will not 
affect the answers and will allow the program to work with smaller structure 
matrices. 

GENERAL NOTES - Multi -harmonic Runs 


As mentioned previously the STARS-2S program has the capability to analyze 
an unsymmetrical loading problem in one submission. In this case the user is 
given a printout option. Upon the insertion of a suitable clue on the data 
card, signifying no intermediate printout, the following printout will be 
provided: 

1. Complete preliminary printout for the first harmonic input (to the 
region joint displacement matrix/see Section 3/)- 

2. Summed stress and deformation data versus the circumferential angle 
0 and meridional distance <J> or s, at the end of program execution. 
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Upon the insertion of a suitable clue on the data card, requesting inter- 
mediate printout^ the following printout will be provided: 

i . 

1. Complete preliminary printout for the first harmonic input (to the 
region joint displacement matrix/see Section 3/)- 

2. Abbreviated stress and deformation data at all segment ends for 
this harmonic. 

3. Segment and region stiffness and load matrices and symmetry checks 
for all other harmonics. 

4. Abbreviated stress and deformation data at all segment ends for all 
other harmonics except the last. 

5. Summed stress and deformation data versus the circumferential angle 
9 and meridional distance <j> or s, at the end of program execution. 


Since the solution of the unsymmetrical loading problem essentially requires 
the consecutive execution of a series of analyses; the program may be in- 
advertently cut off due to an insufficient time estimate. In order to allow 
the user to utilize the good data generated in such an aborted run, a restart 
option has been incorporated into the STARS-2S program. This option is used 
as follows : 

1. After the complete execution of any harmonic, the name of the 
scratch tape presently containing all previously summed harmonic 
data is printed out as scratch tape SAVE l4 or scratch tape SAVE 

, 15. Thus in any computer submission on which a bad time estimate 

can be made, instructions to the operator to save the tapes named 
SAVE 14- and SAVE 15 should be given. These two tapes are the ones 
mounted respectively on logical units l4 and 15. 

2. Upon receipt of an aborted computer run the name of the tape 
containing the last completed harmonic can be obtained from the 
printout. This will be the restart tape; the other can be released. 

3. In order to initiate a restart run the following must be done: 

a) One data card must be changed - the restart clue must be set to 
indicate a restart run. 

b) Operator instructions must be -given to mount the restart tape on 
logical unit 15 ( regardless of the restart tape name ). Another 
scratch tape will be mounted as usual on logical unit l4. 

c) The run will now proceed as normal. In case of a second run 
abort follow the instructions starting with item 2. 
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4. Mote : The restart tape for any given run is not useable more than 

once. Upon the initiation of the computer run this tape becomes a 
scratch tape and additional information (more harmonic data) will 
be written on it. 


GENERAL NOTES - Data debugging . 


The STARS programs have been provided with special separate data debugging 
packages called SATELLITE programs. In order to be able to debug as much 
of a given data deck as possible in one computer submission, the data is 
grouped by inserting special cards, termed "dash-separator cards", appro- 
priately. In order so that additional errors are not made by requiring 
insertion and removal of these cards, the STARS program has been coded to 
accept these dash- separator cards in the input. A dash-separator card is 
shown below: 



As can be seen, a minus symbol is inserted straight across the computer card 
from column 1 through column 80. 

Since the dash- separator cards subdivide the data deck, there exists the pos- 
sibility that a separated data block may be completely omitted (for example 
no kinematic links in a structure). In this case one dash-separator card is 
also omitted . Under no circumstances can there exist two adjacent dash- 
separator cards in a data deck. The SATELLITE programs are described in 


Reference 3- 
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2-3 DETAILED ORDER OF INPUT (See Figures 2-1, 2 - 2 ) 


GENERAL INTRODUCTORY CARDS 


Column Format 


1. Title Card 

A. Alphameric title (submission description) 1-64 l6A4 

2 . Program Control Card 

A. Number of regions to be coupled 1-2 12 

(Max. =19) 

B. Total number of segments 3-5 ■ 13 

(Max.- = 19 x 24 = 456) 

C. Number of Material Property Tables 6-7 12 

(Max. = 10) 

D. Number of Harmonics to be Superimposed 8-9 12 


(Max. =25) 

If the loading is unsymmetric and the 
computer run is comprised of several 
harmonics the following restrictions 
must be followed: 

a) The number of loading conditions 
must be one. 

b) The problem must be linear (LINE 
clue on the master clue card). 

c) If the loading contains thermal 
loads, and if the axisymmetric 
harmonic (n = 0.0) is included, 
then its description must be 
given first. 

If a single harmonic is' to be 
investigated the number of 
harmonics is set to one. 

E. Number of loading conditions for this problem 10-11 12 

The user is able, for each problem submitted, 
to analyze his structure under several inde- 
pendent loading conditions (Max. = 10). The 
number of these loading conditions will 
determine the number of load clue cards 
which will be necessary per segment. If 
the run is to be non-linear or contain 
multi -harmonics to be superimposed, this 
number can only be unxty (l). 
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GRAPHICS CARD 


£ 


0 ANGLES (FOR MULTI- HARMONIC RUNS) 


Repeat for 
each harmonic. 


-T .( 

l 

I 

— £ 


' 0 ANGLE .CONTROL CARD (FOR MULTI- HARMONIC RUNS) 


JOINT LOADS (IF ANY) 


ft 


JOINT LOAD CONTROL CARD (BLANK IF NO JOINT LOADS) 


ft 


ELASTIC SUPPORT STIFFNESSES (IF ANY) 

ELASTIC SUPPORT CONTROL CARD (BLANK IF NO SUPPORTS) 


ft' 


h 


LOAD CLUES, THERMAL 4 DISTRIBUTED LOADS FOR STRUCTURE 


JOINT LOAD CARDS (IF ANY) 


£ 


£ 


[ JOINT LOAD CONTROL CARD (BLANK IF’ NO JOINT LOADS) 

|ELASTIC SUPPORT STIFFNESSES (IF ANY) 

I ELASTIC SUPPORT CONTROL CARD (BLANK IF NO SUPPORTS) 


ft 


W 


BOUNDARY CONDITIONS 


INTER-REGION KINEMATIC LINKS (IF ANY) 


INTER-REGION RINGS (IF ANY) 


JOINT CONTROL DATA CARD 


ft 



£ 


1 


IOTSA-REOIOH KINEMATIC LINKS (IF AM) 


INTRA-REGION RINGS (IF ANY) 


SEGMENT TOPOLOGY 


VALUES (IF NON-LINEAR ) 


TABLE - DISTRIBUTED LOADS (IF ANY) 


TABLE - THERMAL LOADS (IP ANY) 


| LOADING CLUE 


CARD 


ft 


TABLE - WALL CROSS-SECTION GEOMETRY 
| J OR S 
| MASTER CLUE CARD 

j GEOMETRIC DESCRIPTION ' 

MAGIC CONTROL 


SEGMENT IDENTIFICATION CARD 


REGION TOPOLOGY 


REGION IDENTIFICATION CARD 


z : / 

| MATERIAL PROPERTY TABLES 

\ 

| HARMONIC CARDS 

\ 

| PROGRAM CONTROL CARD 

\ \ 

TITLE CARD 

\ 



w. 


\y 


w 


W 


W 


w - 


[V 


iv ■ 


W 




w 


Figure 2-2. Data Sequence 
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GENERAL INTRODUCTORY CARDS (continued) 


Column 


Format 


In the case of a thermal load condition the 
following restrictions will exist: 

a) Only the first . loading condition can 
contain thermal loads. 

b) The output answers will be provided as 

follows: The answers to the first 

mechanical and/or thermal loading con- 
dition will be provided. The answers 
to successive’ mechanical loading con- 
ditions will be superimposed upon the 
thermal load (only) answers of the 
first loading condition. .This is due 
to the fact that the structure stiff - 
'ness matrix generated was that of a 
hot structure with successive applica- 
tions of mechanical loads. 

Note : In a multi-harmonic submission, 

although only one loading condition is 
allowed, all the harmonics may contain 
thermal loadings. 

F. Coupling code 

1. Coupling to occur Code = 1 ; 12 II 

2. No coupling Code = 0 or blank 

Uncoupled runs present a way to size seg- 
ments by use of stiffness- symmetry checks, 
without a full execution of the problem. 

If no coupling occurs program will -give 
only individual stiffness and load matrices 
of each segment; if coupling occurs the 
program will run to completion and give the 
state of stress and deformation of the 
entire structure. In uncoupled runs the 
number of regions is one (l), the total- 
number of segments and the number of seg- 
ments on the first region identification 
card should be the same, and kinematic 
link, discrete ring, boundary condition, 
elastic support, line load, Q angle, and 
graphics cards, and their accompanying 
dash-separator cards should be omitted. 

In addition, the second region introductory 
card (topology), and cards 10 through 12 
of the segment cards are not included in 
an uncoupled run. 
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GENERAL INTRODUCTORY CARDS (continued) 


Column 


Format 


G. Restart clue 13 XI 

1. Not a restart run clue = 0 or blank 

2. Restart run ■ clue = 1 

For a single harmonic the restart clue is 
always equal to zero. For additional in- 
formation on restarting runs see GENERAL 
NOTES - Multi -harmonic Runs (p.2-6). 

H. Intermediate print clue l4 II 

1. Intermediate print requested clue = 0 or blank 

2. No intermediate print clue = 1 

For additional information on intermediate 
print see GENERAL NOTES - Multi -harmonic 
Runs (p. 2 -£), and Section 3- 

I. Graphics clue 1 6 II 

1. Graphic plots requested clue = 1 

2. Graphics not needed clue = 0 or blank 

The graphics option is described later in 
this section with the graphics input card. 

3. Harmonic Cards 

A. The values of the harmonics to be considered in 

a multi-harmonic run (Max. No. = 25). 1-70 14F5.0 

(Up to l4 values on first card, and up to 11 
values on second card if needed.) 

MATERIAL PROPERTY TABLES (Max. = 10 sets) 

As many sets of these cards are used (<10) as there are 
different material property segments in the structure 
to be analyzed. These tables will be used to obtain 
the thermal variation of material properties if thermal 
loadings exist. Thus the range of temperature in this 
table should be greater than that of the thermal loads. 

If no thermal loads exist, the values given in the first 

column of this table will be used, and the rest of the : 

table can be left blank. If there are thermal loads, 

the range of the table is to be considered as that 

between the second and tenth columns . 

1. Identification Card 

A. Material Title (Alphameric) 1-4 A4 

Any name can be made up as long as it is con- 
sistently used on the segment cards to which 
it refers. The same name cannot appear on 
more than one (l) table. 
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MATERIAL PROPERTY TABLES (continued) 


Column 

11-14 


Format 

a4 


B. Type of Table 

One of several possible alphameric clues is 
■written here. These clues serve to size the 
number of cards in the property table, and 
define which properties belong on which card. 

The possible clues are: 

ISOT 
ORTH 
STIF 

Their definitions are provided in Item 2 below. 

2. Material Property Cards 

The material property cards below are given depending 
upon which table type clue is used. If the table 
type clue is "ISOT" (isotropic table): 

A. Temperature values (5 values per card; 2 cards) 5E14.7 

These are the temperatures at which the values 
of material properties will be given. The first 
value in the table must always be the room or 
stress-free temperature, since the material 
properties in only the first column of the table 
will be used in an analysis involving no thermal 
load. The values of temperature in table col- 
umns 2 through 10 must be in algebraically 
increasing order. 

B. Values of Young's Modulus at the given temperatures. 

(5 values per card; 2 cards) 

C. Values of Poisson's Ratio at the given temperatures. 

(5 values per card; 2 cards) 

D. Values of the thermal coefficient of expansion 
at the given temperatures. 

(5 values per card; 2 cards )- 

If the table type clue is "ORTH" (orthotropic table): 

A. Temperature values 5E14.7 

(5 values per card; 2 cards) 

These are the temperatures at which the values of 
material properties will be given. 

B. Values of Young's Modulus in the 6 direction (Eg) 5E14.7 

at the given temperatures. . . 

(5 values per card; 2 cards) 


5El4 . 7 
5El4 . 7 
5El4 . 7 
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MATERIAL PROPERTY TABLES (continued.) 

C. Values of Young's Modulus in the <J> direction (E ) 
at the given temperatures. 

(5 values per card; 2 cards) 

D. Values of the Poisson's Ratio v at the given 

temperatures. ? 

(5 values per card; 2 cards) 

E. Values of the thermal coefficient of expansion in 
the 9 direction (a^) at the given temperatures. 

(5 values per card; 2 cards) 

F. Values of the thermal coefficient of expansion in 
the <j> direction (a ) at the given temperatures. 

(5 values per card; 2 cards) 

G. Values of the Shear Modulus G at the given 

temperatures. ^ 

(5 values per card; 2 cards) 

If the table type is "STIF" (table to be used for 
reinforced shells): 

. A.-G. The values in these locations are the same as 

those above for the "ORTH" clue case, and refer 
to the basic shell. 

H. Values of ring Young's Modulus (E ) at the 
given temperatures. 

(5 values per card; 2 cards) 

I. Values of stringer Young's Modulus (E ) at the 

given temperatures. : 

(5 values per card; 2 cards) 

J. Values of ring thermal coefficient of expansion 
(o^) at the given temperatures. 

(5 values per card; 2 cards) 

K. Values of stringer thermal coefficient of 

expansion (a ) at the given temperatures, 
s 

(5 values per card; 2 cards) 

Note : In a rotated waffle or isogrid construction, 

items H and I, and J and K, refer to the grid 
directions and are respectively identical. 


Column , 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus 

(See General Notes - Data Debugging) 


Format., 
5E14 . 7 

5E14 . 7 

5E14 . 7 

5E14 . 7 

5E14 . 7 

5E1^.7 
5El4 . 7 
5El4 . 7 
5E1^ . 7 

5E14.7 
in 1-80 
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REGION ' INTRODUCTORY CARDS Column Format 

These two cards are placed at the beginning of each region 
data information. Each region contains the following data 
set (see Figure 2-2): a) Two region introductory cards; 

b) data cards for each segment within the region; 

c) ring cards describing the discrete rings within the 
region, if any; and d) kinematic link cards describing 
the kinematic links within the region, if any. 

1. Identification Card 

A. Number of segments within the region (<24 ) 

B. Number of kinematic links between segments 
within the region. 

C. Number of discrete rings between segments 
within the region. 

D. Any alphameric information (region description) 

2. Topology Card (Coupling Orientation) 

A. Region Number 1-5 15 

Number of the region under consideration. 

B. Joint (i) 6-10 15 

Joint associated with i (beginning) end of 
the region. 

C. Joint (j) th 11-15 15 

Joint associated with j (ending) end of 
the region. 

There is no coordinate flow in regions (unless 
1 region = 1 segment), such as that shown for 
segments in Figures 2-3 to 2-9- However, the 
start joint of a region must match with 1 in 
segment numbering, and the end joint must 
match, with the highest segment joint number 
in the region (see Figure 2-12 and page 2-54). 

If 1 region = 1 segment the segment topology 
card will be a dummy card'[l 1 2] (see page 

2-54). 


1-2 

12 

3-k 

12 

5-6 

12 

7-70 

i6a4 
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Figure 2-3. Typical Shell Segment 
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Figure 2-4 a. Ellipsoid 
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+z 



Figure 2-5. Modified Ellipsoid 
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z 



Figure 2 - 7 . General Geometry 
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i in j Specify <}> = 0 



Figure 2-8. Cone 





+z 



Specify r Q 


Figure 2-9- Cylinder 
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SEGMENT CARDS 


Column 


Format 


This sequence of cards is repeated for each segment 
within the region. 

1. Identification Card 


A. 

Segment identification code (see Figures 2-4 to 2-9) 



1 . 

Ellipsoidal or spherical shell 

Code = 11 

1-2 


2. 

Modified ellipse shape 

Code = 12 

1-2 


3. 

Ogival - Toroidal 

Code = 13 

1-2 


4. 

General Geometry 

Code = l4A or 

1-3 



(see Figure 2-7) 

i4b 



5. 

Dummy geometry slot 
(<f> coordinate) 

Code = 15 

1-2 


6. 

Conical - Circular Plate 

Code = 21 

1-2 


The 

plate is treated as a cone with 

zero angle. 



7. 

Cylindrical shell 

Code = 31 

1-2 

B. 

Any 

alphameric information (segment 

description) 

4-67 

"MAGIC" 

Control Card 



A. 

Interval at which final answers are 

to be printed 

1-14 


(in 

radians or inches). 




The 

^-coordinate is defined for all 

geometric 



shapes except the cylinder, cone and plate, 
for which the s coordinate is used. Figures 
2-4 through 2-9 describe these coordinates 
for each shape. 

B. Difference 15-28 

The value recommended depends upon the 
computer used. For eight figure accuracy 
computers it is 1.0 E-6; for the IBM 360 
it is 1.0 E-4. 

C. Integration interval 29-42 

The Runge Kutta numerical integration procedure 
is substantially more accurate than finite 
differencing. An interval of (.03 to. 06) x segment 
size (in radians or inches) should be sufficiently 
accurate for static analyses, (in using a 30 point 
segment table (see p.2-37) there should be at least 
30 integration steps.) 

D. Delta 71-72 

For a fixed-step integration, Delta = 0. 

This card controls the Runge -Kutta numerical 
integration scheme. The suggested values above 
yield accurate results for a fixed-step integration 
method. 


F2.0 

F2.0 

F2.0 

F2.0A1 

F2.0 

F2.0 

F2.0 
16A4 . 

E14.1 


E14.1 


E14.1 


F2.0 
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Calculation of Segment Length 


There is a restriction on the length of the shell segments. Physically, 
the restriction demands that boundary disturbances at one edge be 
distinctly felt at the other edge. This is a consequence of using a 
numerical integration procedure. Since the segment stiffness matrices 
must be symmetric, the calculations involved in obtaining each matrix element 


must be such that a computer round off error never becomes prominent. Limit- 
ing the segment length insures satisfaction of this criterion. This length 
is a function of both geometric shape and segment location within a specific 
geometry. One of the limiting factors is that the ratio of the radii of 
revolution at the initial and final points of a segment be greater than one 

hundredth and less than one hundred. Thus i — I <100 where: 

100 \ r . I 

\ oj / 



This requires smaller segments than will normally be predicted by formula in 
the area of an apex. In addition, note that (cj> = 0) is not an acceptable 
input point (except for the torus-ogive or offset ellipsoid). 


For a cylinder, the segment length parameter, 

A = (l + y) 2 0 As 

should be held to about 4.0. In this expression, "7" is a non-linear param- 
eter. For homogeneous shells: 


■ [* - - ) 


It is zero for a linear problem. 

"f3" is a measure of the rate of decay of a disturbance in the shell. 
"As" is the meridional length. 
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The values of 3 and As for various shell geometries are given below: 



for v ■ 0.3, A s U 


3.11 [ r 0 2 (Uh 2 . 6ht » 3t 2 ! 


2 ' * 
3 - ni [ j|rr] [ <h i * V* * 12 W (b i » h o ♦ *>] ’ 


Approximate formulas can be obtained for near cylindrical regions of generally 
curved surfaces.' The length parameter, 1 

A = (1 + A A* 

should be held to about 4.0. In this expression "y" has the same definition 
as in the cylinder case. 

"A" is a measure of the rate of decay of a disturbance in the shell. 

As 

"A*" = — is the angle intercepted by a meridional arc length "As". 

1 

4 , : • ;■ ■ 

The values of A and As for various shell geometries are given below: 


for v ■ 6.3, A s 4: 


Homogeneous Construction I .‘ 3(1 _ v 2 ) — 


Sandwich Construction - 
Equal Face Sheets 


r s 2 (kh 2 » 6 ht . 3t 2 ) 


i. s 3.11 [ r^IWi 2 ♦ 6bt ♦ 3t 2 )] 


Sandwich Construction - 
Unequal Face Sheets 


. a* 1 - '‘‘V 

(hj ♦ h o ) i * ♦ ♦ h o ♦ t) L P 2 . J 


r r -2 r 

[bffs;] [ li i * h o>‘ ‘ ‘’J 
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_3 

The minimum allowable segment- length is 1 x 10 
(inches or radians according to segment sizing). 


Column Format 

3. Geometric Description Card 

A. Ellipsoid and sphere (Figure 2-4) 



1 . 

Semi-axis perpendicular to Z-direction (a) 

1-14 

E14.1 


2 . 

Ratio of semi- axis in the Z-direction (b) 

15-28 

E14.1 



to (a), 3 = (|) 




3. 

C = offset distance (±) (C = 0.0 if no offset) 

29-42 

E14.1 

B. 

Modified ellipse shape (Figure 2-5) 




1 . 

Axis ratio coefficient (n) 

l-l4 

E14.1 


2 . 

Semi-axis perpendicular to Z-direction (a) 

15-28 

E14.1 

C. 

. .Ogive (Figure 2-6) •' • 




1 . 

r = radius 

l-l4 

ei4.i 


2 . 

C = offset distance (±) 

15-28 

ei4.i 

D. 

General Geometry (Figure 2-7) 




1 . 

Number of pairs of Z versus r Q points (<l4) 

1-2 

12 


2 . 

Z versus r^ points, in pairs, starting 

3-72 

F10.0 



with Z value (7 values per card, including 
first card, for up to 4 cards total). 

1-70 ° r 

F10.0 



In the input table the first Z value is 
taken as Z n = 0, and furthermore Z. >0 

1 l 

(i = 2 -l4) (see Figure 2-7). 



E. 

Cone (Figure 2-8) 




1 . 

Angle 4> in radians (for flat plate, $ = 0). 
Keep in mind that this <p is a constant for 
a given cone and should not be confused 
with the 4> on the MAGIC Control Card. 

l-l4 

E14.1 

F. 

Cylinder (Figure 2-9) 




1 . 

Radius 

l-l4 

E14.1 


4, Master Clue Card 

This card contains a series of clues which de- 
termine the program and table locations to be 
used for the .segment being described. For. a- 
master flow chart of clues and options in the 
. program see Figure 2-1. 
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Column 


Format 


Note: In a multi-harmonic submission the master clue 

card is fixed for all harmonics. Thus for each 
segment. the thermal load clue for each harmonic must 
correspond to the thermal clue on the master clue card. 

A. Material Table Type Clue 1-4 

This clue defines the type of material property 
table to be expected for the segment. This, as 
well as the following clue determines the mate- 
rial properties that will be used in the struc- 
tural analysis for the segment. Thus these- two 
clues should match the two clues used on the 
identification card of the corresponding mate- 
rial property table. As mentioned before on 
page 2-19, the three possibilities are: 

ISOT 

ORTH 

STIF 

B. Material Title 11-14 

This name should be the same as the name which 
appears on the material property table which 
contains the properties to be utilized for 
this segment. 

C. Sheet Clue 21-24 

This clue informs the program as to what kind 
of shell wall crossection to expect. If the 
shell is of single sheet construction, the clue 
to be used is: SING. If the shell wall is an 

equal-size face sheet sandwich, the clue to be 
used is.: EQUA . If the shell wall is a sand- 

wich but the face sheets are not equal, the 
clue to be used is: UNEQ. Finally, if the 
shell segment is reinforced by rings, stringers, a 
waffle, or an isogrid, the clue to be used 
is: BLAN . 

D. Reinforcement Clue 31-34 

This clue describes the type of reinforcement 
that is present on the shell. If the shell is 
purely of single sheet or equal or unequal- size 
face sheet honeycomb construction (no reinforcing), - 
the clue to be used is: THIC. If the reinforce- 

ment consists of rings or stringers or both, 
located along the coordinate axes (0 and <j> or s), 
three clues are possible depending upon the 
basic shell wall construction. If the basic 


A4 


A4 


a4 


a4 
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Column 


Format 


wall construction is a single sheet, the clue to 
be used is Sill . If the basic shell wall is an 
equal-size face sheet sandwich, the clue to be 
used is ST12 . If the basic shell wall is a sand- 
wich but the face sheets are unequal, the clue to 
be used is ST13- If the reinforcement consists 
of a waffle which is rotated at an arbitrary angle 
P from the meridional axis the following three 
clues are possible depending upon the basic shell 
wall construction. If the basic wall construction 
is a single sheet, the clue to be used is RWA1 . 

If the basic shell wall is an equal-size face 
sheet sandwich, the clue to be used is RWA2 . If 
the basic shell wall is a sandwich but the face 
sheets are unequal,- the clue to be used is RWA3 . 

If the reinforcement consists of an isogrid con- 
struction of general angle p from the meridional 
axis (normally p = 30°) the following three clues 
are possible depending upon the basic shell wall 
construction. If the basic wall construction is 
a single sheet, the clue to be used is ISG1. If 
the basic shell wall is an equal-size face sheet 
sandwich, the clue to be used is ISG2 . If the 
basic shell wall is a sandwich but the face sheets 
are unequal, the clue to be used is ISG3 - Two 
other clues are available, namely ST10 and RWAF, 
if the user wishes to input his own stiffness 
constants. These constants may represent any 
wall construction as long as. the basic Hooke's 
Laws used with the clues are appropriate to 
describe the construction to be considered. The 
Hooke's Laws used with these clues are given 
under the description of segment card set 6. 

Note: The reinforcement described in the seg- 

ment cards is closely spaced reinforcement which 
will be smeared over the segment. Discrete rings 
at segment ends in a region are described at the 
end of all segment data for that region. 

E. Thermal Clue 4-1-44 A4 

This clue describes the type of thermal problem 
which exists in the segment. The user is re- 
minded that if there . is a thermal loading on 
the structure, subsequent mechanical load cases 
will be superimposed on the hot structure (see 
page 2-17). If there is no thermal load on the 
segment, the clue to, be used is NOTH . If the 
thermal loading on the segment is of a general, 
standard type, that is if there is variation 
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Column 


Format 


of temperature through the thickness as well as 
in the coordinate directions, the clue to he, 
used is THST . If the thermal load is such that 
the variation is all in the coordinate directions, 
and there is no thermal variation through the 
thickness, the clue to he used is THCN . The last 
clue concerns a shell which is inhomogeneous in 
the meridional direction. This is not really a 
thermal problem at all, but merely a manipulation 
of the material property tables. If a structure 
has a wide variation in material properties in the. 
meridional direction, without this last option one 
must take short segments of constant properties 
for analysis. With this option, however, the 
property variation is placed in the material 
property table, and expressed on the segment as 
a function of temperature. No thermal loads are 
calculated, however, and the temperatures are 
only used to interpolate for material properties 
as integration is progressing along the segment. 

Thus continual variation of properties in the 
meridional direction is accommodated. The clue 
for this option is THIN. 

Note : The program cannot accommodate material 

properties varying in the hoop direction. Thus 
if, in a multi -harmonic submission, a thermal 
load must be described by several harmonics, 
the n = 0.0 information must be provided first. 

The hot shell material properties are then ob- 
tained for this case, and used for all other 
harmonics. 

F. Stress -free Temperature 51 -60 

The value of the temperature (usually room tem- 
perature) at which the segment has no thermal 
stresses or distortions induced, is provided 
here. This is the temperature at which the 
shell was manufactured. If the analysis consists 
of multi-harmonic thermal loads and this value is 
input as other than zero (as would be the case if 

if ATs are input as loads ) , all the harmonic 
amplitude Ts should be of such magnitude that this 
temperature can be subtracted from each to obtain 
the AT in each harmonic. If there is to be no 
thermal analysis, this value is not used and can ' 
be set to zero (O.O). 


E10.1 
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Column 


Format 


G. Non-Linear Clue 

If the analysis is to be linear the clue is LINE. 
For a non-linear analysis see Reference 1 Section 
7 and use the clue NPHI . This latter clue is only 
applicable in a single harmonic (n = 0.0) run. 


A4 


H. Table control - Number of points in each of the 71-72 12 

following tables. 

This can vary from 2 to 30 depending upon the 
shell geometry and loading. For a linearly 
varying geometry and/or loading only 2 input 
points would be required. These two points 
would be the end. points. For more general 
loading and/or geometry a larger number of 
points are required. In particular, each 
abrupt change is specified by two points. 

One should use as many points as necessary 
(up to 30) in order to completely describe 
the problem, rather than using very short 
segments. 

5* Table of 4> or s Values 

A. Initial, intermediate and final values of 5E14.7 

<p or s. Each point requires l4 columns on a 
card and thus there can be 5 values per card 
and up to 6 cards to make a total of up to 
30 points. The points to be specified are 
the beginning point of the segment, any point 
of discontinuity, and the end point of the 
segment. The input must be consistent with 
item H of the previous card. 

6. Table of Wall Crossection Geometry 

The contents of these cards (up to 6 cards per 
item below) are dependent upon the clues reg- 
istered' on the Master Clue Card. If the shell 
to be described contains no reinforcing, the 
pertinent clue is item 4-C, the Sheet Clue. The 
geometry is input and the stiffnesses are cal- 
culated by the program (see Figure 2-10). The 
input is presented below as a function of the 
Sheet Clue. 

If the Sheet Clue is SING (single sheet con- 
struction): 

A. Initial, intermediate and final values of wall 5E14.7 

thickness (h. ) at points defined by table of <}> 
or s values. 
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Figure 2-10 Shell Section Properties 




Column 


Format 


If the Sheet Clue is EQ.UA (equal-size face sheet 
sandwich) : 


A. Initial, intermediate and final values of face 
sheet thickness (h. ) at points defined by table 
of or s values. 

5El4 . 7 

B. Initial, intermediate and final values of core 
thickness (t) at points defined by table of <j> 
or s values. 

If the Sheet Clue is UNEQ (unequal-size face sheet 
sandwich) : 

5E14.7 

A. Initial, intermediate and final values of inner 
face sheet thickness (h.) at points defined by 
table of | or s values. 

5E14.7 

B. Initial, intermediate and final values of core 
thickness (t) at points defined by table of <p 
or s values. 

5E14.7 

C. Initial, intermediate and final values of outer 
face sheet thickness (h ) at points defined by 
table of <f> or s values. 

5E14.7 


If the shell is reinforced, the Sheet Clue will be 
BLAN. In this case it is the following, or Reinforce- 
ment Clue (item 4 d) which will determine the contents 
of card series 6. For the reinforcement cases the 
geometry can be complex and varied, since all types 
of reinforcing are to be included. The reinforced 
shell input is presented below as a function of the 
Reinforcement Clue. 

If the Reinforcement Clue is ST11 (single sheet 
reinforced by rings and/or stringers): 


A. Initial, intermediate and final values of the 5E14.7 

torsional stiffness in the <p direction (GJ ) at 

points defined by table of <j> or s values. 9 

B. Initial., intermediate and final values of the 5E1A.7 

torsional stiffness in the 6 direction (GJg') at 

points defined by table of <p or s values. 

C. Initial, intermediate and final values of 5E14.7 

stringer area (A ) at points defined by table of 

<j> or s values. ? 

D. Initial, intermediate and final values of ring 5E14.7 


area (A^) at points defined by table of f or s 
values. 

E- Initial, intermediate and final values of 

stringer eccentricity (measured inwards from 
base shell centroid as positive) at points 
defined by table of ^ or s values. 
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' Column 

F. Initial, intermediate and final values of ring 
eccentricity (measured inwards from base shell 
centroid as positive) at points defined by 
table of <f> or s values. 

G. Initial, intermediate and final values of 
stringer moment of inertia (about base shell 
centroidal axis) at points defined by table 
of (f> or s values. 

H. Initial, intermediate and final values of ring 
moment of inertia (about base shell centroidal 
axis) at points defined by table of <js or s 
values. 

I. Initial, intermediate and final values of 
stringer spacing at points defined by table 

of 4> or s values. (Do not set to zero if no stringers. ) 

J. Initial, intermediate and final values of ring 
spacing at points defined by table of | or s 
values. (Do not set to zero if no rings.) 

K. Initial, intermediate and final values of base 
shell wall thickness (h. ) at points defined by 
table of 4> cr s values. 

If the Reinforcement Clue is ST12 (equal face sheet 

sandwich reinforced by rings and/or stringers): 

A. through J. The items contained on these cards are 10 sets of 
those described for the ST11 clue above. 

K. Initial, intermediate and final values of base 
shell face sheet thickness (h. ) at points de- 
fined by table of ^ or s values. 

L. Initial, intermediate and final values of base 
shell core thickness (t) at points defined by 
table of $ or s values. 

If the Reinforcement Clue is ST13 (unequal face 

sheet sandwich reinforced by rings and/or stringers): 

A. through J. The items contained on these cards are 10 sets of 
those described for the ST11 clue above. 

K. Initial, intermediate and final values of base 
shell inner face sheet thickness (h. ) at points 
defined by table of <f> or s values. 

L. Initial, intermediate and final values of base 
shell . core thickness (t) at points defined by 
table of <j> or s values. 
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Format 
5El4 . 7 

,5El4 . 7 

5El4 . 7 

5El4 . 7 
5E14 . 7 
5El4 . 7 

5E14 . 7 
5El4 . 7 

5E14 . 7 

5El4 . 7 
5El4 . 7 

5El4 . 7 



Column 


Format 


M. Initial, intermediate and final values of base 
shell outer face sheet thickness (h ) at points 
defined by table of | or s values. 

If the Reinforcement Clue is RWA1 (single sheet 

reinforced by a waffle rotated at an arbitrary angle 

from the meridional direction): 

A. Initial, intermediate and final values of waffle 
grid area at points defined by table of I or s 
values. 

B. Initial, intermediate and final values of waffle 
grid eccentricity (measured inwards from base 
shell centroid as positive) at points defined 
by table of (j) or s values. 

C. Initial, intermediate and final values of waffle 
grid moment of inertia (about base shell centroidal 
axis) at points defined by table of <j> or s values. 

D. Initial, intermediate and final values of waffle 
grid spacing at points defined by table of <j> or 
s values. 

E. Initial, intermediate and final values of waffle 
grid rotation angle, 0 , (in radians from the 
meridional direction) at points -defined by table 
of | or s values. 

F. Initial, intermediate and final values of base 
shell wall thickness (h. ) at points defined by 
table of I or s values. 


" 5E14.7 

5El4 . 7 
5El4 . 7 

5El4 . 7 
5El4 . 7 
5E14 . 7 

5El4 . 7 


If the Reinforcing Clue is RWA2 (equal face sheet 
sandwich reinforced by a waffle rotated at an 
arbitrary angle from the meridional direction): 

A. through E. The items contained on these cards are 5 sets of 5E14.7 
those described for the RWA1 clue above. 


F. and G. The items contained on these cards are 2 sets of 5E14.7 

those described for the ST12 clue above as items 
K. and L. 

If the Reinforcing Clue is RWA3 (unequal face sheet 
sandwich reinforced by a waffle rotated at an 
arbitrary angle from the meridional direction): 

A. through E. The items contained on these cards are 5 sets of 5E14.7 
those described for the RWA1 clue above. 


F. through H. The items contained on these cards are 3 sets of 5E14.7 
those described for the ST13 clue above as items 
K. through M. 
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Format 


If the Reinforcing Clue is ISG1 (single sheet 
reinforced by a general angle isogrid construction): 

A. through D. The items contained on these, cards 
are identical to those described for the RWA1 
clue above, but with reference to the isogrid. 

E. Initial, intermediate and final values of the 
isogrid angle, p, (in radians from the 
meridional direction, see Reference 6 
Appendix A). For the normal isogrid, the 
angle is 30° (input in radians). 

F. . Initial, intermediate and final values of base 

shell wall thickness (h. ) at points defined by 
table of 4> or s values. 1 

If the Reinforcing Clue is ISG2 .(equal face sheet 
sandwich reinforced by a general angle isogrid 
construction): 

A. through E. The items contained on these cards 5 sets of 5E14.7 

are those described for the ISG1 clue above. 

F. and G. The items contained on these cards are 2 sets of 5E14.7 

those described for the ST12 clue above as 
items K. and L. 

If the Reinforcing Clue is. ISG3 (unequal face sheet 
sandwich reinforced by a general angle isogrid 
construction): 

A. through E. The items contained on these cards 5 sets, of 5E14.7 

are those described for the ISG1 clue above. 

F. through H. The items contained on these cards 3 sets of 5E14.7 

ane those described for the ST13 clue above 
as items K. through M. 

If the Reinforcing Clue is ST10 the following Hooke's 
Laws will be used by the program for. the description 
of the shell wall: 

N e = K -n e 0 + K i2 e <t>, ' c n- k e " N T0 
o T o 


4 sets of 5El4 . 7 
5El4 . 7 

5El4 . 7 



2-4-2 


Column 


Format 


M 0 D ll k 0 + D 12 k 4> + C ll e 6 " ^0 

o 

% = r, 22 k (|) + D 12 k e + C 22 £ <p " 

o 


M , 


= -M 


e<p 


“ 2D 33 k 4>0 


Therefore the input is (see Ref. 6 Appendix A): 


A. 

Initial, intermediate and final values of K^ 
at points defined by table of <f> or s values! 

5El4 . 7 

B. 

Initial, intermediate and final values of 
at points defined by table of <f> or s values. 

5E14 . 7 

C. 

Initial, intermediate and final values of 
at points defined by table of <}> or s values. 

5E14 . 7 

D. 

Initial, intermediate and final values of 
at points defined by table of <f> or s values: 

5E14 . 7 

E. 

Initial, intermediate and final values of 
at points defined by table of <f or s values. 
(Should be input as negative for sign convention.) 

5El4 . 7 

F. 

Initial, intermediate and final values of 
at points defined by table of f or s values. 
(Should be input as negative for sign convention.) 

5E14.7 

G. 

Initial, intermediate and final values of 
at points defined by table of <f> or s values? 
(Should be input as negative for sign convention. ) 

5E14.7 

H. 

Initial, intermediate and final values of 
at points defined by table of | or s values: 

5El4 . 7 

I. 

Initial, intermediate and final values of C^ 
at points defined by table of 4> or s values! 

. 5El4 . 7 

J. 

Initial, intermediate and final values of C ^ 
at points defined by table of <J> or s values! 

5El4 . 7 


If the Reinforcing Clue is RWAF the following Hooke's 
Laws will be used by the program for the description 
of the shell wall: 


N e “ K n e 9 + K i2 e <j> " c ii k 9 c i5 k <*> n t0 

o o 

N 4> = K 22 e <j> + K 12 E e " c i5 k 0 " C 22 k $ ' N T4> 
T T o o 


N 4 >e K 33 Y <^e o " 2 c i6 
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Column 


Format 


M e 

M <f> = D 22 k <j> + D 12 k e + c i5 e 6' + C 22 e 4> " ^4, 

o .0 

M <f>e = -2D 33 k 4>e + c i 6 Y <f>e o 

Therefore the input is (see Ref. 6 Appendix A):. 

A. through J. The items contained on these cards 
those described for the ST10 clue above. 

Initial, intermediate and final values of C.. at' 
points defined by table of <p or s values. ^ 

Initial, intermediate and final. values of at 
points defined by table of 4 or s values. 

7- Loading Clue Card (For first harmonic only in a 
multi-harmonic submission) 

The contents of this card are numerical clues which 
alert the program to the types of loads that exist 
on the segment. In addition) if the clue indicates 
that some load does not exist, the appropriate cards 
in series 8 which would ordinarily contain the 
numerical values of this load are omitted from the 
sequence. - 

The series of cards 7 and 8 are repeated for the 
number of loading conditions indicated on the Pro- 
gram Control Card (item E) up to a maximum of 10. 

If for one of these loading conditions, no load 
exists on the segment, then a blank Loading Clue ; 
Card is inserted in the sequence, at this point 
and the load cards are entirely omitted. For 
instance if there are four load patterns to be 
investigated, and for the second pattern some 
segment is completely unloaded, the card sequence 
for that segment would be: 

Loading Clue Card for pattern one (card 7) 

Load Values for pattern one (card 8) 

Blank Card 

Loading Clue Card for pattern three (card 7) 

Load Values for pattern three (card 8) 

Loading Clue Card for pattern four (card 7) 

Load Values for pattern four (card 8) 


K. 

L. 


i 

are . 10 sets of 


- D k e + D k + C e Q + C £ 

y o o 




5El4 . 7 
5E1L . 7 
5El4 . 7 


Column Format 


uiThe -appropriate clues are as follows: 

A. Thermal Clue 

1 II 

If there are no thermal loads (Item is 
NOTH) the clue number is. zero (0). 

If there is a standard thermal variation through 
the thickness (Item HE is THST ) the clue number 
i s four ( 1+ ) . 

If the temperature is constant through the 
thickness (item He is THCN) or if the in- 
homogeneous option is used, (item UE is THIN) 
the clue number is one (l). 

. . If a thermal loading does exist on the structure, 

then the stiffnesses matrix is thermal dependent. , 

In this case successive mechanical load conditions 
. .'.Vi will all be superimposed on the hot structure as 
discussed previously. 

B. Circumferential Load Clue ( f a ) 2 II 

0 

If there are no circumferential loads , then 
the clue number is zero (0). - 

If there are circumferential loads, then the 
clue number is one (l). 

C. Meridional Load Clue (f, ) 3 II 

<P 

If there are no meridional loads , then the 
clue number is zero (O). 

If there are meridional loads , then the clue 
number is one (l). 

D. Normal Load Clue (f^j *4 II 

If there are no normal loads, then the clue 
number is zero (0). 

If there are normal loads , then- the clue 
number is one (l). 

E. Circumferential Moment Load Clue (m Q ) 5 II 

0 

If there are no circumferential moment loads, 
then the clue number is zero (0). 

If there are circumferential moment loads, then 
the clue number is one (l). 
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Column 


Format 


F. Meridional Moment Load Clue (m, ) 6 II 

<p 

If there are no meridional moment loads , then 
the clue number is zero (0). 

If there are meridional moment loads , then the 
clue number is one (l). 

G. Any alphameric information (load description) 7-70 l6A4 

8. Table of Applied Loads (see Figures 2 -11a, b for 
sign convention). Loads are for the. initial har- 
monic only if the submission is multi -harmonic. 

The appropriate card sequence is given below as a 
function of the Loading Clues on card 7* If the Thermal 
Clue is one (l): 

A. Initial, intermediate and final values of the 5E14.7 

temperature of the shell at points defined by 
table of iji or s values. (These values will be 
used either for a thermal problem where there 
is no thermal variation through the thickness 
{Clue = THCN}, or to calculate varying material 
properties along the shell for an inhomogeneous 
problem {Clue = THIN}.) 

If the Thermal Clue is four (4): 

A. Initial, intermediate and final values of the 5E1U.7 

temperature T-^ at points defined by table of 
4> or s values. (The subscripts "nm" indicate 
temperature location - see below. ) 
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Format 


B. Initial, intermediate and final values of the 
temperature T. at points defined by. table of 
<t> or s values . 

C. Initial, intermediate and final values of the 
temperature T at points defined by table of 
i(i or s values . 

D. Initial, intermediate and final values of the 
temperature T at points defined by table of 
4 > or s values . 

If the Thermal Clue is zero (0), the above cards are 
omitted . 

If the Circumferential Load Clue is one ( 1 ) : 

E. Initial, intermediate and final values of the 
circumferential loads f at points defined by 
table of <f> or s values. 

For a discussion of distributed loads see 
Reference k, Appendix A. 

If the Circumferential Load Clue is zero (0), cards E 
are omitted. 


If the Meridional Load Clue is one (l): 


F. Initial, intermediate and final values of the 
meridional loads f at points defined by table 
of <(> or s values . ^ 


If the Meridional Load Clue is zero (0) , cards F are 
omitted. 


If the Normal Load Clue is one (l): 

G. Initial, intermediate and final values of the 
normal loads f at points defined by table of 
4 > or s values . ^ 

If the Normal Load Clue is zero (0), cards G are omitted. 

If the Circumferential Moment Load Clue is one (l): 

H. Initial, intermediate and final values of the 
circumferential moment loads m^ at points 
defined by table of 4 > or s values. 


5E14.7 


5E14.7 


5E14.7 


5E14.7 


5E1U.7 


5El4 .7 


5El4 .7 
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Column 


Format 


If the Circumferential Moment Load Clue is zero (0), 
cards H are omitted. 

If the Meridional Moment Load Clue is one (l): 

I. Initial, intermediate and final values of the 5E1H.7 

meridional moment loads m^ at points defined 
by table of ij or s values; 

If the Meridional Moment Load Clue is zero (0), cards 
I are omitted. 

9- Table of assumed Meridional Membrane Force for 
axisymmetric, single harmonic, non-linear problems. 

If item tG on . the Master Clue Card is NPHI : 

A. Initial intermediate and final values of N^ 5E11.7 

at points defined by table of I or s values. 

If Item 4 g on the Master Clue Card is LINE, card 
sequence 9 is omitted. If the problem under con- 
sideration is non-linear, then the stiffness 
matrix depends upon N(j> and only one load condition 
may be run. In addition the problem may not be 
multi -harmonic. (See Program Control Card.) 

Non-Linear Analysis 

The present STARS-2S shell program is capable of considering non-linear 
effects, but under restricted conditions. The only non-linear behavior that 
can be currently analyzed is that due to axisymmetric loading, and the only 
mode of deformation that is considered under this loading is the axisymmetric 
mode. The iteration technique which is not automated, can be accidentally 
made to diverge. To ensure convergence, care must be taken to satisfy the 
segment sizing parameters. A non-linear analysis would then use the 
following steps: 

l) The first approximation to N is obtained from a linear solution, 
which may yield larger results than the actual non-linear final 
value. 
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2) After the preliminary values are obtained, the structure segments 
should be resized in the areas where local non-linear behavior is 
suspected, using the non-linear sizing parameter (7) described on 
pages 2-31 to 2-33. 

3) If non-linear effects are locally large, the value for ¥ obtained 
from a linear solution (step l) may dominate the loading terms in 
a non-linear analysis, and cause oscillations. To eliminate this 
possibility, the first ¥ approximation should be lower than the 
value predicted by a linear analysis (90 to 50 percent depending 
on the suspected magnitude of non-linearity). If oscillations are 
encountered at any stage, they can be eliminated, and' a trend 
toward convergence reestablished by using smaller' values of N. 

4) With this step the procedure is repeated. A non-linear analysis 
is made and the N output is compared with the barred (assumed) 
quantities. This procedure is carried out until convergence is 
reached. For the magnitude of non-linearities considered by the 
present theory about 2 iterations are usually sufficient. 

Column Format 

10. Stress Clue Card 

This card contains a series of clues which are used to 
identify the proper Hooke's Law to be used for stress 
calculations. If the Reinforcing Clue (item 4 d) is 
THIC, ST10, or RWAF, the sets of cards 10 and 11 are 
omitted . 

A. Meridional Stress Inner Edge Clue 1-4 A4 

This clue informs the program as to what kind ... . y 

of construction exists at the meridional inner 
edge of the shell segment. If the construc- 
tion is part of the basic shell, the clue is: 

SHEL . If the construction is part of the 
rotated waffle reinforcement, the clue is: 

WAFF . If the construction is part of a 
stringer, the clue is: STRI . If the con- 

struction is part of the isogrid, the clue 
is: ISGR . 

B. Meridional Stress Outer Edge Clue 11 -l4 A4 

This clue informs the program as to what kind - ... 

of construction exists at the meridional outer 
edge of the shell segment. The . same possibi- 
lities as in Item A exist, and the possible 
clues again are: 

SHEL . 

WAFF 

STRI 

ISGR 
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' Format 


C. Hoop Stress Inner Edge Clue 

This clue informs the program as to what kind 
of construction exists at the hoop inner edge 
of the shell segment. If the construction is 
part of the basic shell, the clue is: SHEL. 

If the construction is part of the rotated 
waffle reinforcement, the clue is: WAFF . 

If the construction is part of a ring, the 
clue is: RING. If the construction is part 

of . the isogrid, the clue is: ISGR . 

D. Hoop Stress Outer Edge Clue 

This clue informs the program as to what kind 
of construction exists at the hoop outer edge 
of the shell segment. The same possibilities 
as in Item C exist* and the possible clues 
again are : 

SHEL 

WAFF 

RING 

ISGR 

11. Reinforced Stress Calculation Table 

The contents of these cards are dependent upon the 
Reinforcing Clue (item 4 d). If this clue is THIC, 
ST10, or R.WAF, the sets of cards 10 and 11 are 
omitted. 

If the Reinforcing Clue is RWA1, RWA2, RWA3, ISG1, 
ISG2, or ISG3: 

A. Initial, intermediate and final values of 
extreme inward distance to reinforcement 
edge at points defined by ^ or s table. This 
is a signed (±) value, measured from the base 
shell centroid to the extreme point of shell 

, or reinforcement, positive inwards. 

B. Initial, intermediate and final values of 
extreme outward distance to reinforcement 
edge at points defined by ^ or s table. This 
is a signed (±) value, measured from the base 
shell centroid to the extreme point of shell 
or reinforcement, positive inwards. 


21-24 A4 


31-34 A4 


5El4 . 7 


5El4 . 7 
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( 


If the Reinforcing Clue is ST11, ST-12,' or ST13: 

A. Initial, intermediate and final values of 
extreme inward distance on 0 side to reinforcement 
edge at points defined by 4) or s table; 

This is a signed (±) value, measured' from 
the base shell centroid to the extreme 
point of shell or ring, positive inwards. _ 

B. Initial, intermediate and final values of 
extreme outward distance on 0 side to reinforcement 
edge at points defined by or s table. 

This is a signed (±) value, measured from 
the base shell centroid to the extreme 
point of shell or ring, positive inwards. 

C. Initial, intermediate and final values of 
extreme inward distance on cp side, to reinforcement., 
edge, at points defined by 4 or s table. 

This is a signed (±) value, measured from 
the base shell centroid to the extreme 
point of shell or stringer , positive 
inwards. 

D. Initial, intermediate and final values of 
extreme outward distance- on cp side to reinforcement 
edge at points defined.. by 4 or s table;. 

This is a signed (±) value, measured from 
the base shell centroid to the extreme 
point of shell or stringer , positive 
inwards. 


12. Segment Topology Cards 

A. Segment number 

Number of the segment under consideration. 

B. Joint (i) 

Joint associated with i^* 1 end of the segment 
(TIC). 

C. Joint ( j ) 

Joint associated with the j end of the 
segment (STOP). 


Since within a region the segments are all 
singly connected, the segment joint numbers 
should be in adjacent numerical pairs. That 
is, if joint (j) is 6, joint (i) could only 
be 5 or 7. This is true only within a region. 


1-5 


6-10 


11-15 


5KL4.7' 


5E14.7 


5El4 . 7 


5El4 . 7 


15 

15 


15 
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In addition, the initial joint of each 
region must be 1 in segment topology 
numbering, and the final joint of each 
region must be the last (highest) number 
in the segment topology numbering (see 
Figure 2-12) . The coordinate $ or s 
increases from TIC to STOP, i to j . 

The user is again advised to see 
Figures 2-3 to 2-9. 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in 1-80 

(See General Notes - Data Debugging) 


INTRA-REGION DISCRETE RING CARDS 


These cards, if any exist (region introductory card, 
item 1C), are placed at the end of all the segment 
data for the region. They contain the following 
information for each ring (in groups of 3 cards): 


1. Ring Location and Property Card 

A. Segment joint number to which the ring is 1-2 12 

attached. If there is also a radial 
discontinuity at this location necessitating 
a kinematic link, the joint number of the 
ring is the independent joint number of the 
link (see Figure 2-13a). 


B. 

Ring 

extensional stiffness (EA) 

3-i6 

El4. 7 

C. 

Ring 

axis 

bending stiffness about centroidal y 
(Ely). See Figure 2-l4. 

17-30 

El4 . 7 

D. 

Ring 

axes 

cross-bending stiffness about centroidal 

(El ). 
xy 

31-44 

E14.7 

E. 

Ring 

torsional stiffness (GJ) 

45-58 

El4 . 7 

F. 

Ring bending stiffness about centroidal 

x axis (El ). 

' x 

59-72 

E14.7 


2. Ring Geometry Card (see Figure 2-l4) 


A. Ring thermal coefficient of expansion (a^). 


Ring centroidal radius (r c ). 


E12. 5 


B. 


1-12 

13-24 


E12.5 






z 


z 

4 


u- 

No Links Necessary- 




Independent joint 
of kinematic link 
and therefore Ring joint 
number . 


Kinematic Link Necessary 


(a) Segment or Simple Region Rings 



4 





Independent joint for all 
kinematic links and 
therefore Ring joint 
number. 


(b) Possible Region Multi-connected Rings 


Figure 2-13. Discrete Ring Topology 
2-57 
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Format 


C. Distance (±) between ring centroid and shear 25-36 

center in the x direction (x ). 

c 

D. Distance (±) between ring centroid and shear 37-48 

center in the y direction (y ). 

E. Offset (±) of ring shear center from attached 49-60 

shell joint (see item LA.) in the x direction (x). 

F. Offset (±) of ring shear center from attached 61-72 

shell joint (see item 1A) in the y direction (y) . 

Note ; The shear center (S ) is the origin of the 
x, y axes shown in Figure 2-l4. 

Ring Thermal Description Card (See Figure 2-l4) 

A. Distance (±) between ring shear center and l-l4 

the closest (innermost) radially measured 

ring point x-j. (in x = radial direction). 

B. Distance (±) between ring shear center and the 15-28 
farthest (outermost) radially measured ring 

point, Xq (in x = radial direction). 

Note : Items A & B cannot both be input as zero. 

C. Innermost point temperature (Tj). 29-42 

Note : For a harmonic n ^ 0.0 this is the 

amplitude of a cosine expansion Tj cos n6. 

D. Outermost point temperature (Tq). 43-56 

Note : For a harmonic n^O.O this is the 

amplitude of a cosine expansion Tq cos n0. 

E. Ring stress-free temperature (room or manufac- 57-70 
turing temperature (see Note for Master Clue 

Card item F) . 

The number of sets of discrete ring cards 
must be equal the number specified in item 
1C, of the region introductory card. 

Note : Thermal loads on the rings are input 

on the above cards. Mechanical loads on 
the ring must be input as a load distribution 
on the segments connecting to the ring. A 
higher intensity loading may be distributed 
near the segment ends connected to the ring. 


E12.5 

E12.5 

E12.5 

E12.5 

E14.7 

E14.7 

El4. 7 
E14.7 
E14.7 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D 
(See General Notes - Data Debugging) 


minus in 1-80 



Column 


Format 


INTRA-REGION KINEMATIC LINK CARDS 

These cards, if any exist (region introductory card, 
item IB), are placed at the end of all the discrete 
ring data for the region. They contain .the following 
information: 

A. Joint ( j ) dependent joint . , 1-2 12 

B. Joint (i) independent joint ■ > 3-^ 12 

For intra-region kinematic ■ links these joints 
must be in consecutive dece'nding order. That 
is, joint ( j ) should always 1 be greater than' 1 
joint (i) by one. 

C. Angle y in radians (see Figure 2-15) 5-19 El^.7 

7 cannot equal 0 or ir. . 

The angle y describes the orientation of the 
link; it is the inclination angle of the link 
from the vertical (Z axis). 

The number of kinematic link cards must 1 
equal the number specified in item IB, 
of the region introductory card.' 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus In 1-80 

(See General Notes - Data Debugging) 

REGION JOINT CONTROL DATA 

These cards are placed at the end of all the data for 
all regions. 

1. Joint Control Data Car'd 1 i 


A. 

Number of region joints 

1-5 

15 


Total number of region joints (Max. = 20). 



B. 

Number of discrete rings 

6-10 

15 


Total number of discrete rings '.between * 

regions. 



C. 

Number of kinematic links' 

11-15 

15 


Total number of kinematic links between 
regions. 
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Column 

Discrete Ring Cards (inter-region, if any exist) ' 

A. Ring Location and- Property Card ' 

a. Region joint number to which the ring is 1-2 

attached. If there are also radial 
discontinuities at this location neces- 
sitating kinematic links, or multi- 
connections, the joint number of the ring . i 

is the independent joint number of all • 
the links (see Figure 2-13b). , / 

b. Ring extensional stiffness (EA) 3-l6 

c. Ring bending stiffness about centroidal 17-30 

y axis (Ely). See Figure 2-l4. . . 

d. Ring cross-bending stiffness about ,< 31-44 

centroidal axes (El ). 

' xy 

e. Ring torsional stiffness (GJ) ■ •> 45-58 


f. Ring bending stiffness about centroidal x 
axis (El ). 


x 


B. Ring Geometry Card (see Figure 2-l4^..' 

a. Ring thermal coefficient of expansion 

( V- i 

r' 

b. Ring centroidal radius (r ). , 

° / 

c. Distance (±) between ring centroid and 
shear center in the x direction (x ). 

/ 

d. Distance (±) between ring /centroid and 

shear center inthe y direction (y ). 

/ ^ 


e« Offset (±) of ring shear/ center from 

attached shell joint (see item Aa) in the 
x direction Cx). 


f. Offset (±) of ring shear center from . 
shell joint (see item/Aa) in the y 
direction (y). j> 

Note : The shear center (S) is the origin 

of the x,y axes shown in Figure 2-l4. 


59-72 


1-12 

13-24 

25-36 

37-48 

49-60 

61-72 


Format 

12 

E14.7 

E14.7 

E14.7 

E14.7 

E14.7 

E12.5' 

E12 . 5' 
E12.5 

E12.5 

E12.5 

E12.5 



C. Ring Thermal Description Card (See Figure 2-lk) 


Column 


Format 


a. Distance (±) between ring shear center 1-lk 

and the closest (innermost) radially 

measured ring point, Xj (in, x = radial 
direction). 

b. Distance (±) between ring shear center 15-28 

and the farthest (outermost) radially 

measured ring point, Xq (in x = radial 
direction) . 

Rote: Items a & b cannot both be input as zero. 


c. Innermost point temperature (Tj). 29-k2 

Rote : For a harmonic n ^ 0.0 this is the 
amplitude of a cosine expansion Tj 

cos n6. '• 

b 

d. Outermost \point temperature (Tq) . k3-56 

Rote : For a harmonic n ^ 0.0 this is the 
amplitude of a cosine expansion Tq 

cos n9. 

e. Ring stress -free temperature (room or 57-70 

manufacturing temperature (see Rote for 

Master Card item F) . ' 

The number of sets of discrete ring cards 
must equal the number specified in item 
IB, of the Joint' Control Data Card. 

' .4 

Rote : Thermal loads on the ring are input 
on the above cards v . Mechanical loads on 
the ring must be input as line loads or 
moments on the shell joint (item Aa) 
assigned to the ring. 

D-A-S-H S-E-P-A-R-A\T-0-R C^A-R-D 
(See General Rotes - Data Debugging) 

\ 

3. Kinematic Link Cards (inter -region, if any exist) 

A. Joint (j) dependent joint V 1-2 

B. Joint (i) independent joint \ 3-^ 

For kinematic links between regions there 
are no restrictions upon joint numbering. 


minus in I 


Elk. 7 
Elk. 7 

Elk. 7 
Elk. 7 
Elk. 7 


-80 

12 

12 
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Column 


Format 


The only restriction is that between successive 
kinematic link data cards the (j) joint entry 
should be in increasing order (not necessarily 
consecutive). For example: 


j_ 

Must be ( 3 

in 1 4 

increasing j 9 

order. ( 12 

C. Angle 7 in radians (see Figure 2-15) 

7 cannot equal 0 or it. 



5-19 


The angle 7 describes the orientation of the 
link; it is the inclination angle of the 
link from the vertical (Z axis). The number 
of kinematic link cards must equal the number 
specified in item 1C of the joint Control 
Data Card. 


E14.7 


D-A-S-H S-E-P-A-E-A-T-O-R C-A-R-D minus in 1-80 

(See General Notes - Data Debugging) 

Boundary Condition Cards (Joint Data - one card 
per region joint) 


A. 

Joint Number 

1-2 

12 

B. 

Joint component conditions on: 




1) Ap 

3-4 

F2.0 


2) or ^ (see Figure 2-l6) 

5-6 

F2.0 


3) ^ ^ 

7-8 

F2.0 


■P- 

CD 

9-10 

F2.0 

C. 

Angle CL in radians 

11-24 

E14.1 


To be used only in conjunction with a 2 or 
3 code. 

Note : There must be as many boundary con- 

dition cards as there are joints as indicated 
in item 1A of the Joint Control Data Card. 
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Figure 2-15. Kinematic Link 
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There are 4 different codes that are used 
to prescribe joint component conditions. 
They are:. 

a. 0 = no displacement allowed. 

b. 1 = displacement allowed in the 

indicated direction. 


A^. and A are rotated through an 
angle or (tt/2 - a) and become A^ 
and A, respectively, while a 
displacement is allowed in the A^ 
direction. iv 

A^ and A^are rotated through an angle 
or (tt/2 - a) and become A^ and A 
respectively, while a displacement 


is allowed in the 


direction. 


See Figure 2-l6 for a geometric 
explanation of codes 2 and 3- 

When using rotation codes: 

Code 2 can exist only as A^ coding. 

Code 3 can exist only as A^ coding. 

Codes 0 and 1 can appear in either 
column 4 or column 6, in addition 
to columns 8 and 10. Thus, there 
are twelve possible boundary con- 
ditions when rotation codes are 
used . (a = (fi for table below.) 



Ap, n e free 


Ap, fixed 


A^, fixed 
free 


Ap free 
fixed 


Free edge (possible 
to apply shear and/ 
or membrane loads) 

Aq = 0, normal 
support (possible 
to apply membrane 
load) 

An = 0, membrane 
support (possible to 
apply shear load) 

' 1,2, 3,1 

1,2, 0,1 

1,0, 3,1 

0,2, 3,0 

0,2, 0,0 

0,0, 3,0 

0,2,3,1 

0,2, 0,1 

0,0, 3,1 


1 , 2 , 3,0 


1 , 2 , 0,0 


1 , 0 , 3,0 
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Apex boundary conditions: 


Since the closed apex (<f> = 0) angle 
is not acceptable input, the apex 
boundary conditions must be simulated 
at a small <p angle. These boundary 
conditions vary per Fourier harmonic, 
and are as follows: 


n 



a r 

n e 

- angle 

0.0 

0 

0 

3 

0 

' a = <p 

1.0 

1 

2 

0 

l 

a = <p 

>2.0 

0 

0 

0 

0 - 

— 


For a further discussion- see 

Reference 6, Appendix C. , . 

Note : In a multi-harmonic submission for a 

shell containing an apex, some approximation 
must be made since boundary conditions . are 
not changed. If there is no load near the 
apex, and the actual membrane stress resultants 
are unknown, it can be left totally free. 

General Notes : 

1) To establish a datum for measuring 
displacement, free body motion must be 
eliminated from the structure. This should 
be accomplished by suitably applied boundary' 
conditions. 

2) The ability of a dependent joint in a 
kinematic link to prescribe motion independ- 
ently should be removed by setting all 
boundary conditions of that joint to zero. 

See pages 2- 62 and 2- 78. 

3) Any joint to which an elastic' support ' 
is attached must be left free to move 
(elastically) in that direction. Therefore 
the boundary code for that direction is 1. 

k) In an axisymmetric analysis, if no 
torsion loads are applied, a faster solution 
will result if all torsion degrees of 
freedom are removed. 
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Column 


Format 


D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D 
(See General Notes - Data Debugging) minus in l-=80 


ELASTIC SUPPORT DATA (for initial harmonic only in a 
multi -harmonic submission) 

1. Support Control Data Card 

A. Number of Elastic Supports < 20 1-4 14- 

Note: Elastic supports can be applied 

only to region joints. 

B. Any alphameric information (support 5-68 16 a4 

description) 

Note : If there are no elastic sup- 

ports for the structure for. the initial 
harmonic, card 1 of the SUPPORT DATA 
is blank and card set 2 is omitted. 

2. Support Stiffness Cards (as many as in 
item 1A above) 


A. 

Supported Region Joint Number 


1-5 

15 

B. 

Stiffness applied to region 
torsional direction (T,,). 

joint 

in 

6-19 

El4 .7 

C. 

Stiffness applied to region 
Z direction (T,.^). 

joint 

in 

20-33 

E14.7 

D. 

Stiffness applied to region 
radial direction (T^)* 

joint 

in 

34-47 

El4 .7 

E. 

Stiffness applied to region 
moment direction (T^). 

joint 

in 

48-6l 

El4 .7 


Notes : l) The use of these support 

stiffnesses is varied. Since they are 
merely added into the structure stiffness 
matrix, they can represent actual elastic 
(spring) supports at region joints. If 
the analyst anticipates the existence of 
a ring at any region support, but has no 
design details and is thus unable to 
provide the ring information necessary 
on the ring cards, a preliminary estimate 
of the ring may be input here by stiff- 
nesses. This option can also be used 
to simulate additional attached structure. 
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Column 


Format 


2) The boundary condition codes for 
the supported .joints should be 1 for 'the 
supported directions. 

•D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus 

(See General Notes - Data Debugging) 

JOINT LOAD DATA (for initial harmonic only in a multi- 
harmonic submission) 

1. Load Control Data Card 

A. Number of Joint Loads 1-4 

Total number of joint loads in analysis. 

(Line loads can-., only be applied to region 
joints. ) 

B. Any alphameric information (load description) 5-68 

Note : If there are no Joint Loads for the 

structure for the initial harmonic, card 1 of 
the JOINT LOAD DATA is blank and card set 2 
is omitted. 

2. Joint Load Cards (as many as in item 1A above) 

A. Load Condition Number 1-5 

Number of loading condition in which the line 
load exists. (See Program Control Card item 
2E. For a multi -harmonic submission this 
number is always unity. ) 

B. Row Identification 6-10 

The identification is the location of the 
degree of freedom at which the load is 
applied. This is obtained by counting the 
non-zero codes entered in the Boundary 
Conditions Cards, starting with Joint 1; 

T, Z, R, Og, Joint 2 ; T, Z, •••••• etc . y 

and stopping at the joint and degree of 
freedom where the line load is to be 
applied. The location number of this 
degree of freedom is the information 
necessary. 


in 1-80 


14 


16a4 


15 


15 
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Column 


Format 


C. Applied Joint Load 11-24 El4 .7 

The input is,2irr times the running load 
in lb. /in. In tfie particular case of the 
axial axisymmetric load, this is simply 
the net force. For sign convention see 
Figures 2-11, 2-l6 and 2-17. 

D-A-S-H S-E-P-A-R-A-T-O-B C-A-R-D minus in 1-80 

(See General Notes - Data Debugging) 

MULTI-HARMONIC DATA (This data set (items 1-4) is 
omitted for single harmonic submissions. Note 
that only one loading condition is allowed for 
multi -harmonic submissions.) 

This data set consists of distributed segment load 
data, elastic support data, and line load data for 
the structure, per harmonic. 

1. Segment Distributed Load Data (This data is 
provided in order of the processed segments 
for the whole structure .) 

A. Loading Clue Card for next harmonic (see 
previous discussion in segment card set 7) 

a. Load clues for temperature, circum- 1-6 6ll 

ferential, meridional, normal, circum- 
ferential moment, and meridional moment 

loads, for next harmonic. 

Note: Variations in the types of 

loadings as well as magnitudes are 
allowed per harmonic. The thermal 
clue however is fixed and must remain 
constant for all harmonics . Thus if 
the first harmonic clue is 4 for 
example, and if the other harmonics 
are not thermally loaded, the Thermal 
Clue must remain 4, and zero (or stress- 
free) temperature values must be entered 
for other harmonics on the load cards. 

b. Table Control (see segment card 4 item H). 71-72 12 

This value does not change per 

harmonic. 
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Figure 2-17L Line Loading for Harmonic n = 1.0 


Column 


Format 
5El4 . 7 


B. Table of Applied Loads for next harmonic 
(see segment card set 8). These cards' .. • 
contain the next harmonic loads as de- 
scribed previously in the segment data, 
(card set 8) according to the Load Clues 
above . 

Note: The above data sets (1A .and IB) are 

repeated until all the segments in the 
structure are completed for the next 
harmonic. 

2. Elastic Support Data for next harmonic 

A. Support Control Data Card for next harmonic 

a. Number of Elastic Supports < 20 

b. Any alphameric information 

Note : If no supports exist in this 

harmonic the above card is blank, 
and card set 2B (below) is omitted. 

B. Support Stiffness Cards for next harmonic, 
(as many as in item Aa above) 

a. Supported Region Joint Number 

b. Stiffnesses applied to region joint in 
torsional, Z, radial, and moment direc- 
tions for next harmonic (see previous 
Elastic Support Data Cards). 

3. Joint Load Data for next harmonic, 

A. Load Control Data Card for next harmonic 

a. Number of Joint Loads 

b. Any alphameric information 

Note : If no joint loads exist in this 

harmonic the above card is blank, and 
card set 3B (below) is omitted. 

B. Joint Load Cards for next harmonic (as 
many as in item Aa above) 

a. Enter unity in column 5 


. 1.-4 
5-68 


1-5 

6-6l 


1-4 f 

5-68 ‘ 


1-5 


l4 

i6a4 


15 

4E14.7 


. 4 

‘ i6a4 


15 
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Column 


Format 


b. Row identification (see previous Joint 6-10 15 

Load Card description) 

Note : Locations of line loads may change 

from harmonic to harmonic. 

c. Applied Joint Load 11-24 El4 . 7 

The input is 2-rrr 0 times the running load 
in lb. /in. for this harmonic (see Figure 
2-17). 

Note: The above data sets (l, 2, 3) are 

repeated for every harmonic until all the 
information is provided for every harmonic 
on the Harmonic Cards (see p. 2-18). 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in 1-80 

(See General Notes - Data Debugging) 

4. Circumferential Angle Data 

A. Theta Angle Control Card 

a. Number of theta angles (l < n < 36). 1-2 12 

The number of theta angles that is 

required for a full description of the 
circumferential variation of the output 
quantities. 

b. Any alphameric information 3-66 l6A4 

B. Theta Angle Cards (<8) 

a. Theta angles in radians (Max. =36). 5E14.7 

The theta angles at which the output is 
to be provided. A maximum of 36 angles 
is allowed making up at most 8 cards. 

GRAPHICS DATA (This data is omitted if the clue on the 
Program Control Card item I was set to zero.) 

1. Graphics Control Card 

This card is divided into two data fields. The 
first field, 1-11, is for representing the item 
to be plotted, while the second field, 21-56, is 
for choosing the theta angle for which the output 
amplitude is to be plotted. 
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Column 


Format 


A. Circumferential Displacement Clue 1 

Clue = 1, plot u 

Clue = 0 or blank, do not plot u 

i 

B. Meridional Displacement Clue 2 

Clue = 1, plot v 

Clue = 0 or blank, do not plot v 

C. Normal Displacement Clue : 3 

Clue = 1, plot w 

Clue ' = 0 or blank, do not plot w - 

D. Transverse <j> Shear Stress Clue 4 

Clue = 1, plot r £<p - 

Clue = 0 or blank, do not plot t , 

£<P 

E. Hoop Stress (inner face) Clue - x . 5 

Clue = 1, plot oq i n 

Clue = 0 or blank, do not plot °q i n 

F. Meridional Stress (inner face) Clue • 6 

Clue = 1, plot a <j) i n 

Clue = 0 or blank, do not plot in 

G. In-plane Shear Stress' (inner face) Clue . 7 

Clue = 1, plot T ( j ) e 'i n 

Clue = 0 or blank, do not plot t • . 

<p 0 m 

H. Transverse 0 Shear Stress Clue 8 

Clue = 1, plot X £0 

Clue = 0 or blank, do not plot i ^ 

I. Hoop Stress (outer face) Clue 9 

Clue = 1, plot o e out 

Clue = O' or blank, do not plot 

J. Meridional Stress (outer face) Clue 10 

Clue =. 1, plot °<j ) ou t 

Clue = 0 or blank, do not plot a, 

’ <p out 

K. In-plane Shear Stress (outer face) 'Clue . 11 

Clue = 1, plot 

Clue = 0 or blank, do not plot qu ^. 

L. - Theta angle location clues (or load case 21-56 

clues - see Note 2 for single harmonic 
submissions) . 

Each item, starting in column 21, .rep- .. 
resents an angle on the Theta Angle Card in 


II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

3611 
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Column 


Format 


the- order of that input. Thus column 21 
represents the first theta angle, column 
22 represents the second theta angle and ■ 
so on up to the number input (Max. = 36 ). 

If the amplitude of the output variables 
of interest (chosen items from A through 
K) are to be plotted for a given angle, 
unity is inserted in that corresponding 
column. If no plots are to be made, the 
corresponding column is to be left "Blank. 

Notes : l) Since in this way up to 396 

plots per segment could be made, an ar- 
bitrary maximum number of plots was set 
at 56 . This number is to be used as a 
maximum on the product of the chosen 
items to be plotted and the number of 
angles at which the plots are to be made. 

2) In a single harmonic sub- 
mission, the field 21-^0 on the Graphics 
Card represents the load cases (up to a 
maximum of 10 ) for which the chosen 
amplitudes (from columns l-ll) are to be 
plotted. Unity is inserted if the load 
case is to be plotted, and a zero or 
blank if not. The field 31-56 is left 
blank. The upper limit of .56 plots per 
segment still holds. 

D-A-S-H S-E-P-A-R-A-T-O-R C-A-R-D minus in 1-80 

(See General Notes - Data Debugging) t 

NEW PROBLEM CARDS 

In summary, in order to repeat a problem' or to do several problems, 
the following options are available. 

1. Loading changes - The present program allows up to ten separate 
loading cases on single harmonic submissions." This applies to 
distributed and line loads. No backup, cards are needed for this 
purpose. (See input information.) 

2. Harmonic changes - The program allows for the superposition of up to 
25 harmonics. There can only be one load case with this option. 

The option is valid for coupled or uncoupled runs. No -backup cards 
are needed for this purpose. (See input information.) 

3. - Full data duplication - The number of full data decks that can be 

backed up on one submission is -limited only by time requirements. 
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It is hoped that the user is now able to use the STARS-2S program to good 
advantage. It is a powerful tool., which will increase in value to the user 
as he uses, it.- One of the more complex areas of usage is the description 
of topology, especially when involved with rotation codes and joint loads. 
An illustrative example of a Y joint representation is therefore presented 
below. (See Figures 2-l8 and 2-19 for the structure and idealization.) 

The idealized structure contains four regions and two kinematic links. The 
joints are numbered from 1 to 7- Membrane loading is applied to joints 1 
and 5 and the structure is supported by membrane action at joint 7* All 
regions must be coupled and we are interested in the 4^ harmonic only.' 


The second card in each region description (topology card) is as follows: 


Region 

1 

2 

3 

4 


Joint (i) 

Joint (j) 

I 

1 

2 

2 

3 

4 

5 

6 

7 


The Joint Control Data- card would contain a 7 in column 5 and a 2 in column 15. 

In this example, the restraints at joints 1, 5 and 7 must be rotated from the 
fixed (global system) to a local system such that membrane action may be 
applied. 'In addition, joints 3 and 4, and 6 and 3 are to be coupled with 
kinematic links. Thus, the motions of joints 4 and 6 are dependent upon the 
motion of joint 3- This dependence will be insured by using 2 kinematic link 
cards and setting the displacements of joints 4 and 6 equal to zero. It should 
be noted in this particular case that the motion of joints 4 and 6 is not being 
equated to zero, but rather, the ability to prescribe motion independently is 
being removed. The required data has the following appearance. 


Kinematic 

Link Cards (2 cards) 


Joint (j) 

4 

6 

Joint 

3 

3 

ii) 

Angle 

7 + 7T 

7 

i 

I Note: In a double, link of the 

l type shown in Figure 2-19, one 
l joint must be consistently 
' independent (joint 3 in example). 
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Figure 2-l8. Y- Joint (Distorted Geometry) 
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B. Boundary Condition Cards (7 cards) 

Joint _T_ _Z_ R_ _0_ Angle 

1 1 2 ' 3 1 

2 1111 

3 11.11 

4 0 0 0 0 ; 

5 1 2 3 1 

6 oooo 

7 0 0 3 l” a 3 

The Load Control Data card would contain a 2 in column 4. 

The external membrane loads (one load problem assumed) are applied to the 
structure through the Joint Load' Cards which, in this example, would appear 
as (2 cards): 

Row Load • 

2 2TT r o x Membrane Load (lb. /in.) 

l4 ^ 7rr 0 x Membrane Load (lb. /in.) 


Problem No. 
1 
1 
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SECTION 3 


OUTPUT INFORMATION 


'The output of the STARS -2S program is straightforward, however a 
description is in order since the user should learn the significance 
of the various checks that are provided. It is important to point 
out that the output of the program will include a print-out of the 
input data. This gives the user the opportunity to check whether 
or not the input data was correct. In the detailed description of 
the complete output which follows, the user should refer to the output 
of the problems in Section 4 as an example. 

The title page of the output contains all the data from the General 
Introductory Cards, prominently placed, and needs no comment. The 
next page of the output contains the first region Identification 
Card in the center. The following output is then presented for each 
segment in this region (in order of appearance): 

1. Contents of segment Identification Card. 

2. Contents of MAGIC Control Card 

3. Contents of Geometric Description Card (Cards) 

4. Contents of Master Clue Card 

5. The material property table used for the segment 

6. Crossection description table 

7. Temperature load table (if any) 

8. Distributed load tables (printed per loading condition) 

9. Contents of Non-Linear Cards (if any) 

10. Segment influence coefficients (MAGIC output) 

11. Segment stiffness matrix 

12. Stiffness matrix symmetry check 

13. Segment load matrices 

14. Radius of revolution at i"^ and j*'* 1 ends of segment 
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Item 12, the stiffness matrix symmetry check, is a check upon, the 
validity of segment sizing and the accumulation of round-off ■ error. 

For perfect symmetry to exist, it is necessary to have zeros above 
the main diagonal, and zeros or ones below the diagonal. The amount 
of error induced by improper sizing or round-off is related to the 
amount that the off-diagonal terms in the lower triangle differ 
from unity. An attempt should be made to keep the- upper limit on 
this difference at one percent (maximum number in lower triangle 
should be 0.1010... E Ol). 

As mentioned previously, items l-l4 are repeated for all the segments 
within region one. The radius of revolution at every joint should be 
checked at corresponding joints of adjacent segments to make sure 
that proper coupling has been specified. At this stage in the output, 
the topology of the segments within the region, and the description 
of the intra-region rings and kinematic links is presented. 

Next the region matrices are presented. Given in order are: the': 

region stiffness matrix, the stiffness matrix symmetry check, and 
the region load matrices. Again the numerical round-off, evident 
in the symmetry check, should be kept to a maximum of one' percent 
(0.1010... E 01 in lower triangle). The output to this point, that 
is, sets of items l-l4, segment topology and links, and the region 
matrices, are now repeated as a group for each region within the 
structure. When this is completed, the region, topology is presented. 

The next items to be provided by the output are the descriptions of 
the inter-region rings and kinematic links and the boundary conditions. 

At external points of the structure these are physical boundary 
conditions. At internal points they merely state the fact that no 
restraint exists and the joint in question is free to move. The last 
column in this set, gives the angle a, which is zero unless a rotation 
code is indicated . It is important to refer to Figure 2-6 once more and 
point out that a represents a rotation of the coordinate system. The 
final item in this output section presents the elastic support description. 
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There are a variety of errors that can be made in submitting input data. 
The STARS-2S program is set so that as an error is detected, 

(input or in the matrix calculations) the program is stopped. 

Therefore, to avoid delays in getting an answer, the SATELLITE-1S 
program was created to debug input data, and to check for other 
possible inconsistencies. The use of this program is described in 
Reference 3, and all the STARS error messages are contained in the 
listings therein. 

After having corrected any errors, so that it is now possible for the 
program to run to. completion, the problem solution can be discussed. 

The first item provided which represents the solution of the problem 
is the structure flexibility matrix. This matrix should also be 
checked for symmetry. After the reduced flexibility matrix, the 
applied line loads and the region end deflections are presented 
per problem. The deflection values are given for the region joints 
in numerical order, starting from joint 1. 

From this point the output will depend upon the intermediate' print 
clue. If this clue is set for full output the following items will 
be printed. First there will be a reproduction of some input data 
previously discussed. Next, there appears a diagram which outlines 
the format of the results tabulated per harmonic, at the beginning 
and end points of each segment. This tabulation is explained in 
Table 3-1- Below this tabulation will appear the numerical results 
at the segment ends. The following special cases may arise: 

1. The 1,2 element will be zero for cylinders or cones. 

2. The 2,3 and 2,4 elements will be equal for linear analyses. 

All the above results will be for harmonic amplitudes until the 
final harmonic is completed. They will be separated by the segment 
and region harmonic stiffness and load matrices for each harmonic. 

At that time a different and more extensive tabulation diagram will be 
presented as shown in Table 3-2. The numerical results will then be 
printed for the prescribed print interval, rather than just segment ends. 
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TABLE 

3-1. 

Tabulated Intermediate Output Information 

PHI (RAD. OR IN.) 
U 
V 
W 

j 

cp or s at which calculations are made 

circumferential displacement 

' 

meridonal displacement 
normal displacement 

DEGREES 
Q PHI 
Q THETA 
OMEGA THETA 
OMEGA PHI 


cp expressed in degrees 

Q , meridional transverse shear stress resultant 
cp' 

Qq, circumferential transverse shear stress resultant, 
j circumferential rotational displacement (rad) 
meridional rotational displacement (rad) ■ 

PRINT INERVAL 
J PHI 
N THETA 
M THETA 


interval at which final answers are printed 

J , effective meridional transverse shear stress resultant 
9 

N , circumferential force resultant 

M 0 , circumferential bending moment resultant 

STEP 

J PHI STAR 
N PHI 
M PHI 


numerical integration step size 

J*. non-linear effective tranverse shear stress resultant 
, 9 

N , meridional force resultant 
9 

M , meridional bending moment resultant 

R ZERO 
T PHI THETA 
N PHI THETA 
M PHI THETA 


. r , radius of revolution 

o 

Tp 0> effective in-place shear stress resultant 
N^ 0 , in-plane shear stress resultant 
M^gj twisting moment resultant 

BASE THICKNESS 
N TEMPERATURE THETA 
N TEMPERATURE PHI 
M TEMPERATURE THETA 
M TEMPERATURE PHI 

basic shell wall thickness 

N^g, circumferential temperature force resultant- load 
Nrp , meridional temperature force resultant load 
M^,g, circumferential temperature moment resultant load 
M^, meridional temperature moment resultant load 


3 - 








TABLE 3-2. Tabulated Final Output Information 


1.1 PHI (RAD. OR IN.) 

2 . 1 EPSILON THETA 

3.1 U 

4,1' V 

5.1 W 

6.1 OMEGA THETA 

7.1 OMEGA PHI 

1.2 DEGREES 

2.2 EPSILON PHI 

3.2 Q PHI 

4.2 J PHI 

5.2 Q THETA 

6.2 TAU ZETA PHI = Q/t 

7.2 TAU ZETA THETA = Q/t 

1.3 PRINT INTERVAL 

2.3 GAMMA PHI THETA 

3.3 K PHI THETA 

4.3 N THETA 

5.3 M THETA 

6.3 SIGMA THETA IN 

7.3 SIGMA THETA OUT 

1.4 STEP 

2.4 K PHI 

3.4 J PHI STAR 

4.4 N PHI 

5.4 M PHI 

6.4 SIGMA PHI IN 

7.4 ' SIGMA PHI OUT 


9 or s at which calculations are made 
e Q , circumferential strain 
circumferential displacement 
meridional displacement 
normal displacement 

U£j, ' circumferential rotational displacement (rad) 
meridional rotational displacement (rad) 


9, expressed in degrees 
e , meridional strain 

Q , meridional transverse shear stress resultant 

9 

J , effective meridional transverse shear stress 
4 resultant 


Qq, circumferential transverse shear stress resultant 
T C9 = meridional transverse shear stress 

T £9 = %/t' c i rcum; f eren 4ial transverse shear stress 


interval at which final answers are printed 

\p9> shear strain 

kqjg, specific twist 

Nq, circumferential force resultant 

Mg circumferential bending moment resultant 

, 

a„. , circumferential stress on inside fiber 
win 

a 9out’ stress on outside fiber 


numerical integration step size 
meridional curvature 

J* , non-linear effective transverse shear stress 
4* resultant 

N , meridional force resultant 

9 

M__, meridional bending moment resultant 
, meridional stress on inside fiber 

9m 

<3 , , meridional stress on outside fiber 

90 ut’ 
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TABLE 3-2. (Continued) 

1.5 

2.5 

3.5 
^,5 

5.5 

6.5 

7.5 

R ZERO r , radius of revolution 

o’ 

K THETA kg, circumferential curvature 

T PHI THETA T^g, effective in -plane shear stress resultant 

N PHI THETA N^g, i n -Pl a ne shear stress resultant 

M PHI THETA M^g, twisting moment resultant 

TAU PHI THETA IN T in-plane shear stress on inside fiber 

TAU PHI THETA OUT T ^ 0 OU £, in-plane shear stress on outside fiber 

1,6 

2,6 

3.6 
M 

5.6 

6.6 
7,6 

. 

BASE THICKNESS basic shell wall thickness 

N TEMPERATURE THETA N^q, c f rcum ferential temperature force resultant load 

N TEMPERATURE PHI ^Tcp’ mer 13-i ona l temperature force resultant load 

M TEMPERATURE THETA M^q, circumferential temperature moment resultant load 

M TEMPERATURE PHI M„ , rtier f^f° na f temperature moment resultant load 

SIGMA F IN °Fin’ ^uber-von Mises-Henoky effective stress 

SIGMA F OUT CT Fout’ Hut) ® r " von Mises-Hencky effective stress 
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and in a Fourier sum form per theta angle. The following new special 
cases may arise: 

1. The 6,2 and 7,2 elements will not be calculated for non- 
sandwich construction and will be printed as zero. 

2. The 6,6 and 7,6 elements will not be calculated for 

-- ■■■■ • reinforced constructions and will- be printed as zero. 

3. No stresses will be calculated for the ST10 or RWAF options. 

4. The 6,3 and 7, 3 and the '6,4 and 7,4 elements will 
actually be in the waffle grid directions for RWA1, RWA2, 
or RWA3 options (see Appendix B). 

; 5. The 6,5 and 7,5 elements will not be calculated at 

- reinforcement edges (they will be calculated in the basic 
•w . shell) and will be printed as zero. . . 


If the intermediate print is turned off, all the print discussed above 
pertaining to harmonic amplitudes will not be available. A restarted 
run will always provide the initial print items 1-14, and will then 
continue from the restarted, harmonic depending upon the intermediate 
print clue.. This completes the program output. The user is reminded 
.to .check continuity, of stress resultants across segment boundaries 
where applicable. 


In utilizing the program, it is frequently necessary to relate applied 
edge loads to the net forces across a section. The relation between 
forces in the fixed (global) coordinate system and any rotated 
coordinate system is given by(References 2 and 4): 


[ jp \ (n) 

1 

0 

0 

0 

I ^ ( 

0 

— s 4> 

-c4> 

0 

) f r ( 

0 

+c<i> 

—Sp 

0 

\ M / 

0 

0 

0 

1 



(3 


for any harmonic "n" 


(see Figure 2-l6b) . 



The relations between the net resultant external loads and the 
magnitudes of distributed edge, loads are:, £! . ; . 



+x +y 


Antisymmetric Loads, (n = T) , 

Here, the integration of the distributed 
forces is not obvious. 

2?r . 2 tt : 

? X * /(- F T (1)r o M ) sin 6 • f (v l V 5 ) cos 6 

o , t 0 , v. .... . 
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When the loaded edge has standard coordinates, Equations (3-2 or 3-3) 
are used to relate the applied edge loads to the new forces across a 
section. In the axis'ymmetric case, this is straightforward: 

~ ( o) 

contributions to F^ are made only by F^ and there. is a similar 
~ (o) 

relation between and F^ . However, in the antisymmetric case, 
there are four unknowns and only two equations. Thus, additional data 
is required. Often, these loads are applied in a region of assumed 

membrane stress; then F^^ = 0 (cylinder) or F^ 1 ^ = 0 (plate), and 

=0 , since these are transverse shear and bending stress resultants. 
When the load is applied at an angle to the fixed (global) coordinates. 
Equation (3-l) is used to transform Equations (3-2 and 3-3) into a form 

which permits evaluation of F^^ and F^^ (for membrane problems only: 

<t> = a, M = F = 0) . 


It is frequently desirable to be able to calculate net forces at a cut 
section, or a built-in edge directly from output values, in order to 
check equilibrium. The net forces in terms of the stress resultants 
(in local coordinates as they appear in the output) are: 
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where the sign is chosen to correspond with the edge ( i or j) on which 
the applied force is desired. 

Equations (3-^ through 3-7) should be used to check overall shell 
equilibrium for unfamiliar geometries since this is a good check on 
the solution to the problem. 



. SECTION 4 

EXAMPLES OF INPUT AND SOLUTIONS 

This section contains the full input and output of an example test 
problem to the STARS -2S program. The ‘problem was chosen to illustrate 
a variety of different inputs rather than from a physical point of view. 

The test problem is presented in Figure 4-1. The clamped spherical 
cap is to be analyzed under a multi-harmonic loading involving both 
distributed and line loads. In addition, a ring is simulated by 
elastic support data as would be the case in an early version of 
design where the final ring geometry is unknown. The shell segments 
are of various wall cross-sections, exemplifying typical input 
requirements. The intermediate print clue has been set for full output 
and graphic displays of some output variables have been requested. 
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STRUCTURAL PERTINENT DATA 


Region 1 

a) Segment 1: Single Sheet 

b) Segment 2: Single Sheet reinforced- by 
waffle 

c) Joint 2 has a ring simulated by 
elastic support input 

d) Segments 1 and 2 are separated by a 
kinematic link 


Region 2 

a) Segment 1 : Equal face sheet sandwich 
reinforced by isogrid 

b) Segment 2 : Unequal face sheet sandwich 
reinforced by rings and stringers 

c) Joint 4 is fixed 

d) Region joints 2 and 3 are separated by 
a kinematic link 


Figure 4-1 Test Problem 
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SEGMENT NUMBER 1 SEGMENT CODE il FIRST PART Of- REG. NO. 



SEGMENT load matrices 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY -UNITS TO SI UNITS 


The International System of Units (Si) was adopted by the 
Eleventh General Conference on Weights and Measures in i960. Conversion 
factors for the units used in this report are given in the following 
table: 


Physical quantity 

U.S. Customary 
Unit 

■ 

SI Unit 
(**) 

Length 

in. 

0.0254 

meters (m) 

Stress modulus 

ksi 

6.895 X 10 6 

newtons/ me ter 2 
(N/m^) 

Stress resultant 

lbf/ in. 

175.1 

newtons/ meter 
(N/m) 

Temperature change 

°F 

5/9 

Kelvin (K) 


* Multiply value given in U.S. Customary Unit by conversion factor to 
obtain equivalent value in SI Units. 

** Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

giga (G) 

10 9 

mega (M) 

10 6 

kilo (k) 

10 3 

deci (d) 

10" 1 

centi (c) 

10- 2 

milli (m) 

10" 3 
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APPENDIX B 


STRESS CALCULATIONS FOR REINFORCED CASES 


Several items must be noted when a stress analysis is performed upon a 
reinforced construction by use of the present program. Some of these 
items can be treated as rules and some require engineering judgement. 
The rules to remember are as follows: 

1. The hoop and meridional stresses output in a waffle reinforcement 
are actually in the waffle rib directions. 


2. In order to obtain grid stresses for an isogrid construction 
the following tranformation "equations are necessary: 



e 


2 


+ e 
£ 

2 


+ 



cos (tt+2(3) + 


\p9 sin (tt+29) 
2 



(^- 1 ) 


Thus from the STARS-2S strain output the appropriate grid strains 
may be obtained, and thence, by using Hooke's Law, the stresses. 
Note that the strains used in Equation (4-l) should be at the 
desired crossection location thus: 


e = e + C.k 
10 1 


(^- 2 ) 
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3 . For a reinforced section, the in-plane shear stress is calculated 
only at a shell location : (if requested) such as point 2 in 
Figure B-l, rather than points 1 or 4. 

4 . For an ST10 or RWAF option the crossection is described by the 

' user in terms of overall stiffness properties. Thus no stresses 
can be calculated by the program for these cases. 

Items requiring engineering consideration in the stress analysis, are 

the following: 

1. If a shell crossection contains materials of different properties, 
either actually or due to differential thermal loading, there is 
no guarantee that the stress at the extreme point is the most 
critical. The analyst should decide whether a check of stresses 
at each location, where material properties change is necessary. 

2. For a structure such as that shown in Figure B-l, the program will 
calculate direct stresses at points 1 , 4 , 3 and the bottom point 
below 3, and in-plane shear stresses at points 3 and the point 
below 3. (Of course this can be modified by inputing different 
offset; distances and stress clues.) While strain is assumed 
linear from point 1 to point 4 , there is a stress discontinuity 
at point 2. This will occur even if the material properties of 
shell and ribs are the same, since the governing Hooke's Laws 
differ (see Reference j). The analyst must again decide whether 
the stresses at point 2 could be more critical than those at the 
extreme points. 
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Any additional stress calculations that the analyst decides to make on 
the basis of items 1 and .2 above, should be always made using the 
program strain and curvature output, and Hooke's Laws for ring, 
stringer or shell (Ref. 7). This can be done automatically in the 
program by, providing the correct Q distance (instead of the extreme) 
and by setting the stress clues for the appropriate Hooke's Laws 
(segment card data sets 10 and ll). 
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Figure B-l Special Stress Case 
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